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CHAPTER  I 
INTRODUCTION 


The  laser  Doppler  Anemometer  (LDA)  has  become  a  valuable  tool 
for  experimental  investigations  in  fluid  mechanics  since  it  can  measure 
the  direction  and  magnitude  of  the  velocity  accurately  and 
nonintrusively.  Today  there  are  commercial  pointwise  LDA  systems 
available  with  specially-designed  slow-moving  traversing  mechanisms  to 
obtain  complete  flow  pattern  in  large  wind  tunnels.  That  means 
for  each  spatial  location  the  experiment  has  to  be  repeated.  Scanning 
LDA  systems  are  needed  for  two  reasons.  First,  having  almost 
"instantaneous"  velocity  profiles,  one  can  investigate  the  flow 
structure.  Secondly,  scanning  LDA  speeds  up  the  acquisition  of  velocity 
information  in  at  least  one  direction. 

This  thesis  is  concerned  with  the  development  of  the  Scanning 
Laser  Doppler  Anemometer  (SLDA)  and  demonstration  of  its  usefulness  in  a 
turbulent  separated  flow  region. 

1.1  Previous  Work: 

There  have  been  a  few  scanning  LDA  systems  designed  and  used  in 
the  past.  The  ones  known  to  the  author  will  be  discussed  here. 

In  1971  a  scanning  LDA  system  was  described  by  Bendick  (1971)  in 
which  he  used  a  translational ly  oscillating  mirror  as  the  heart  of  his 
design.  This  system  was  used  for  instantaneous  velocity  measurements  in 
steady  and  pulsating  water  flow  in  a  glass  tube  of  6mm  I.  D.  with  scan 
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speed  of  approximately  0.4  m/sec.  The  proper  operation  of  this  design 
at  high  scan  speeds  and  large  scan  ranges  is  doubtful  because  of 
inertial  considerations. 

In  1973,  a  two-color  dual-beam  backscatter  scanning  LDA  system 
was  reported  by  Grant  and  Orloff  (1973).  The  scan  was  accomplished  by 
translating  a  lens  in  the  direction  of  the  optical  axis  which  caused  the 
measurement  volume  to  move  in  that  same  direction.  Scan  ranges  from 
60cm  to  200cm  with  variable  scan  rates  from  less  than  1  cm/sec  to 
maximum  of  1.5  m/sec  reported  to  have  been  achieved.  In  later  reports, 
this  system  was  employed  to  investigate  periodic  flow  around  a 
helicopter  blade  as  it  passed  some  fixed  azmuthal  angle  using  a  mineral 
oil  smoke  generator  for  seeding  purposes.  Because  of  low  data  rates, 
measurements  of  instantaneous  velocity  profiles  were  not  possible,  and 
many  scans  were  needed  to  represent  a  plot  of  mean  velocity  with  a 
counter  type  processor.  There  were  no  attempts  to  represent  turbulent 
velocity  fluctuations,  and  fast  scan  rates  were  presumably  difficult  in 
this  translational  system.  For  more  information  concerning  the 
application  of  this  design,  see  Orloff  and  Biggers  (1974)  and  Orloff 
et  al .  (1975). 

In  1976  a  design  description  of  a  backscattered  scanning  system 
was  reported  by  Rhodes  (1976)  in  which  the  scan  range  of  30cm  and  a  scan 
frequency  of  30  Hz  were  reported.  Measurements  of  flow  velocities  at  16 
discrete  positions  along  the  optical  axis  were  made  using  a  large 
rotating  wheel  containing  16  ports,  15  of  which  held  glass  plates  of 
different  thickness.  The  resident  time  at  each  position  is  2  msec. 
This  scanning  system  was  used  to  investigate  the  flow  behind  the  airfoil 
using  kerosine  smoke  particles  as  seeding  material  with  a  measurement 
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volume  dimension  of  0.2  mm  diameter  and  1.3  cm  length.  A  counter  type 
processor  was  used  to  find  the  mean  velocity  for  each  scan.  For  more 
information,  see  Gartrell  and  Jordan  (1977),  and  Meyers  (1979). 

A  new  optical  system  was  reported  in  1978  by  Nakatani  et  al . 
(1978)  which  is  capable  of  making  true  instantaneous  velocity  profiles. 
In  contrast  to  all  previous  designs,  there  is  no  moving  or  scanning 
device.  Spherical  lenses  are  used  to  increase  the  beam  diameter  of  the 
two  horizontal  incident  laser  beams.  After  passing  through  a  converging 
cylindrical  lens,  these  intersecting  beams  form  a  vertical  measurement 
volume  along  a  straight  line.  The  scattered  light  from  this  vertical 
measurement  volume  is  received  by  a  series  of  optical  fibers  connected 
to  individual  photo-detectors,  each  collecting  data  simultaneously  from 
different  positions  of  the  measurement  volume.  Measurements  of  unsteady 
flow  in  a  branch  tube  of  rectangular  cross-section  (2  mm  x  10  mm)  has 
been  reported  in  their  paper.  The  results  clearly  demonstrate  the 
proper  operation  of  this  design  for  small  tubes,  but  instantaneous 
velocity  profiles  in  a  large  tunnel,  using  this  design,  is  expensive  and 
not  practical. 

Durst  et  al .  (1981)  in  their  scanning  system  used  a  rotationally 
oscillating  mirror  to  scan  the  measurement  volume  perpendicular  to  the 
optical  axis  of  the  transmitting  optics.  The  image  of  the  measurement 
volume  for  all  the  positions  fall  on  one  physical  point  on  the 
photomultiplier  (PM-tube).  The  position  transducer  as  well  as  PM-tube 
signals  were  simultaneously  digitized  by  a  two-channel  transient 
recorder,  which  was  triggered  only  when  a  Doppler  burst  was  present. 
The  mean  and  rms  velocity  profiles  via  this  scanning  system  were 
compared  with  the  pointwise  measurements  and  the  conclusion  was  due 
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to  resonance  difficulties  of  the  scanner  mirror  scan  frequencies  were 
limited  to  about  15  Hz.  As  the  system  is  designed,  the  measurement 
volume  scans  an  arc  instead  of  a  straight  line,  and  for  a  given  scanner 
mirror,  fringe  spacing  cannot  be  varied. 

1.2  Objective  of  the  Thesis 

Although  the  scanning  laser  doppler  anemometer  has  various 
useful  applications,  the  entire  work  here  was  motivated  by  the  fact  that 
in  a  separated  flow  region  instantaneous  velocity  gradients  are  sub¬ 
stantially  different  from  the  mean  velocity  gradients.  Mixing  length 
and  eddy  viscosity  models  fail  in  this  region  (Simpson  et  al . ,  1981). 

Mean  velocity  and  turbulence  profiles  have  already  been  mea¬ 
sured,  but  no  satisfactory  model  for  the  backflow  has  emerged. 
(Simpson,  1983).  If  an  almost  instantaneous  velocity  profile  could  be 
obtained,  then  more  detailed  features  of  the  nature  of  the  backflow 
would  be  revealed. 

With  current  technology,  it  is  not  easy  to  determine  purely 
instantaneous  profiles  unless  a  large  number  of  LDA  systems  are  used 
simultaneously.  Therefore,  the  objectives  of  the  present  work  are  as 
follows : 

1.  Development  of  a  rapid  scanning  LDA  system  capable  of  scanning 
at  least  40  cm  or  longer  at  high  frequencies. 

2.  Demonstration  of  the  proper  operation  of  the  above  design  and 
its  usefulness  in  the  separated  flow  occurring  in  the  SMU  wind 


tunnel . 
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1. 3  _ Out  1  in e  of  Thesis 

Chapter  2  describes  the  geometry  and  dimensions  of  the  wind 
tunnel  as  well  as  the  details  of  the  solid  particle  feeder  employed  for 
flow  seeding  in  the  tunnel.  The  conical  manifold  used  for  uniform  flow 
seeding  is  also  described  in  this  chapter. 

Chapter  3  deals  with  the  design  of  the  new  scanning  LDA.  The 
components  of  the  transmitting  and  receiving  optics  are  explained  in 
conjunction  with  a  discussion  of  a  non-dimensional  parameter 
representing  the  number  of  fringes  crossed  by  a  particle  in  this 
scanning  system.  Analysis  of  the  velocity  biasing  is  also  mentioned. 

Chapter  4  is  devoted  to  the  detailed  description  of  the  data 
acquisition  and  processing  procedures  and  logics.  Position  and  velocity 
uncertainty  analysis  as  well  as  the  methods  used  for  testing  the  entire 
system  are  also  explained  here. 

Chapter  5  gives  the  experimental  results  taken  in  a  separating 
turbulent  boundary  layer  occuring  in  the  SMU  wind  tunnel.  The  first  set 
of  results  is  for  demonstrating  the  proper  operation  of  this  new  design, 
and  the  system  is  used  to  map  out  the  flow  in  the  tunnel  to  further 
reinforce  its  capabilities.  Later  in  this  chapter,  space-time 
correlation  coefficients  taken  for  the  first  time  in  a  separated  flow 
are  shown  and  discussed  with  the  hope  of  acquiring  some  more  insight  to 
the  flow  structure.  The  suggestions  for  further  improvements  and 
conclusions  are  mentioned  at  the  last  part  of  this  chapter. 
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CHAPTER  II 

EXPERIMENTAL  APPARATUS 


2.1  Basic  Wind  Tunnel 

The  mainstream  flow  of  the  blown  open-circuit  wind  tunnel  is 
introduced  into  the  test  section  after  first  passing  through  a  filter, 
blower,  a  fixed-setting  damper,  a  plenum,  a  section  of  honeycomb  to 
remove  the  mean  swirl  of  the  flow,  seven  screens  to  remove  much  of  the 
turbulence  intensity  and  finally  through  a  two-dimensional  4:1 
contraction  ratio  nozzle  further  to  reduce  the  longitudinal  turbulence 
intensity  while  accelerating  the  flow  to  test  speed. 

The  upper  wall  is  adjustable  so  that  the  free-stream  velocity  or 
pressure  gradient  can  be  adjusted.  The  side  walls  are  made  of  float 
glass.  The  test  wall  is  constructed  from  1.9  cm  thick  fin-form  plywood, 
reinforced  every  28  cm  with  7.6  x  3.8  x  .6  cm  cross-section  steel 
channel.  Figure  2.1  is  a  side-view  schematic  diagram  of  the  7.62  m 
long,  .91  m  wide  test  section  of  the  wind  tunnel. 

The  active  boundary-layer  control  system,  which  is  described  by 
Simpson,  Chew,  and  Shivaprasad  (1980),  is  used  to  eliminate  preferential 
separation  of  the  curved-top-wall  boundary  layer.  Highly 
two-dimensional  wall  jets  of  high-velocity  air  are  introduced  at  the 
beginning  of  each  of  the  .20  m  long  sections.  At  the  latter  two 
streamwise  locations  the  oncoming  boundary  layer  is  partially  removed  by 
a  highly  two-dimensional  suction  system. 
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The  inviscid  core  flow  is  uniform  within  0.05%  in  the  spanwise 
direction  and  within  1%  in  the  vertical  direction  with  a  turbulence 
intensity  of  0.1%  at  20.9  m/sec.  Figure  2.2  shows  the  free-stream 
velocity  distribution  obtained  along  the  tunnel  center-line  using  the 
single-wire  probe. 

2.2  Solid  particle  Seeder 

The  choice  of  the  particle  seeder  explained  in  this  section  is 
based  on  a  series  of  exhaustive  experiments  conducted  on  the  entire 
system  including  SMU  wind  tunnel  and  counter  type  processor  for  data 
acquisition.  (See  Chapter  IV) 

Adrian  et  al .  (1975)  showed  that  a  signal  with  at  least  45dB 
signal-to-noise  ratio  (SNR)  was  required  for  results  not  prejudiced  by 
the  counter  processor  settings.  Since  determination  of  SNR  practically, 
if  not  impossible,  is  not  easy  for  all  seeding  agents,  a  more  realistic 
approach  was  used  to  find  the  "proper"  seeding  material  in  this  work. 
Here,  "proper"  means  that  the  particles  are  faithfully  following  the 
flow  and  the  45  dB  SNR  criterion  for  counter  processor  is  fulfilled. 

To  find  the  appropriate  seeding  material  for  this  work,  the 
results  (mean  and  RMS  of  velocity)  for  the  pointwise  velocity 
measurements  taken  with  the  scanning  set  up  (explained  in  Chapter  3),  at 
a  given-fixed  streamwise  position,  were  compared  with  the  known 
established  results  (reference  velocity  profiles)  measured  with 
different  arrangement  at  that  same  location  (Simpson,  1978). 
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2.2.1  Seeding  Techniques  Examined 

The  aerosol  generator  at  SMU  was  based  on  the  design  proposed  by 
Echols  and  Young  (1963)  using  dioctyl  phthalate  (DOP)  as  seeding 
material.  A  series  of  experiments  using  this  aerosol  generator  with 
scanning  LDA  arrangement,  operated  in  pointwise  mode,  did  not  give 
satisfactory  results.  It  was  suspected  that  too  many  sub-micron  DOP 
particles  were  present  and  were  producing  much  signal  noise  with  a 
fringe  spacing  of  3.3  pm. 

A  series  of  detailed  tests  on  this  DOP  aerosol  particle 
generator  were  performed  at  the  University  of  Arkansas  by  the  group 
headed  by  Professor  M.  K.  Mazumder  using  the  SPART  aerosol  analyser 
(Mazumder  et  al .  1979).  These  tests  confirmed  the  suspected  factor 
since  only  30%  of  the  particles  were  ever  greater  than  1  pm.  In  a  chain 
of  experiments  a  technique  proposed  by  Seegmiller  (1982)  which  uses 
Paasche  HS  US  air-brush  with  a  mixture  of  2  liters  of  denatured  meth¬ 
anol,  lOOcc  of  warm  water,  and  lOcc  polylatex  spheres  was  examined. 

Air  brushes  were  used  in  various  streamwise  locations,  just 
upstream  of  the  leading  edge  of  the  test  wall  of  the  tunnel  and  upstream 
of  the  contraction,  and  at  different  injecting  angles  with  respect  to 
the  floor  of  the  tunnel.  No  "satisfactory"  results  were  obtained  with 
the  scanning  LDA,  operated  in  pointwise  mode  having  3.3  pm  fringe 
spacing,  in  comparison  with  the  reference  velocity  profiles. 

Various  powders  were  used  with  the  existing  particle  seeder,  ex¬ 
plained  in  the  next  section;  and,  finally,  it  was  found  that  the  baby 
powder  produced  relatively  'satisfactory'  results.  The  estimated 
density  of  the  baby  powder  is  about  0.8  g/cm3  which  is  almost  650  times 
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that  of  air.  The  particle  diameter  is  also  estimated  to  be  in  the  order 
of  few  microns  (<10  um).  Hjelmfet  and  Mockros  (1966)  showed  that,  for 
the  case  where  the  ratio  of  the  seed  density  to  fluid  is  about  1000; 
particles  with  7  ym  diameter  will  give  frequency  response  up  to  812  Hz. 
This  limiting  frequency  increases  with  decrease  in  particle  diameter  and 
when  density  ratio  approaches  one. 

2.2.2  Details  of  the  Seeder  Design 

The  solid  particle  seeder  built  at  SMU  is  similar  in  design  to 
the  Opto-Elektronische  Instrumente  (Karlsruhe,  Germany)  (OEI)  Model  101 
solid  particle  seeder  which  uses  baby  powder  having  particles  greater 
than  1  ym. 

Air  passes  through  a  Sears  air  line  filter-regulator  Model  No. 
282-160232,  and  into  the  seeder  housing.  It  is  important  that  the  air 
be  dried  to  prevent  additional  agglomeration  of  baby  pc-.d'-  part,;V-s. 

The  schematic  of  the  seeder  is  shown  in  Fig.  2.3..  The  outlet  of 
the  housing  is  made  to  contract  gradually.  The  inside  corners  are 
filled  with  sealing  foam  to  streamline  the  interior  surface  and  thus 
prevent  baby  powder  from  piling  up  in  the  housing.  One  side  of  the 
housing  is  made  of  plexiglas  which  enables  one  to  inspect  the  inside 
during  operation.  The  inside  of  the  box  is  accessible  via  a  sealed  top 
lid.  The  box  pressure  can  be  increased  up  to  1  psig  without  significant 
leakage. 

A  1/4  HP,  1725  RPM,  single  phase  AC  motor  is  used  to  rotate  the 
shaft  inside  the  housing  via  a  1  to  1  pully  mechanism.  Four  15.2  cm 
diameter,  fine  Sears  wire  wheels  (Model  No.  9HT64824)  are  mounted  on  the 
shaft  inside  the  box  which  is  being  supported  by  two  self  aligned 
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bearings  outside  the  housing  and  are  fixed  to  the  particle  seeder  frame. 
The  two  side  holes  of  the  box,  from  which  the  shaft  passes  through,  are 
sealed  by  two  rubber  flat  plates  having  a  hole  size  equal  to  that  of  the 
shaft  O.D.  on  each  of  them  and  are  mounted  with  two  collars  on  each  side 
to  reinforce  the  sealing  due  to  higher  inside  pressure  of  the  housing. 

The  baby  powder  is  pushed  up  inside  the  housing  by  an  aluminum 
piston  of  4.78  cm  diameter,  30.73  cm  length,  and  is  scraped  by  wi.e 
wheels  inside.  A  variable  DC  motor,  to  control  the  seeding  rate,  with 
calibrated  tachometer  and  2  to  1,  90°  speed  reducer  rotates  a  10:1  worm 
gear  system  which  in  turn  translates  a  lead  screw  vertically  upward. 
This  lead  screw  pushes  the  piston  inside  the  cylinder  to  feed  the 
housing  with  the  power  (Figs.  2.4  and  2.5). 

Air  flow  rate  through  particle  seeder  is  measured  via  a  Dwyer 

3 

Rotometer  up  to  4  m  /hr  and  is  directed  towards  the  wire  wheels.  The 
tip  velocity  of  wire  brushes  is  approximately  4  m/sec,  which  is  much 
higher  than  the  maximum  jet  speed  of  0.3  m/sec. 

The  outlet  of  the  housing  is  directed  to  the  conical  manifold 
which  is  described  in  the  next  section.  The  qualitative  response  of  the 
seeder,  using  counter  processor,  is  shown  in  Fig.  2.6.  The  estimated 

3 

average  concentration  of  the  baby  powder  is  33  mg/m  and  it  takes  about 
10  -  15  minutes  to  reach  its  almost  steady  state  condition. 

2.3  Conical  Manifold  for  Particle  Dist ribution 

It  is  important  to  seed  the  flow  in  the  wind  tunnel  (for 
measurement  purposes)  as  uniformly  as  possible  so  as  to  minimize  its 
effects  on  two-dimensional  character  of  the  boundary  layer.  The  powder 


Fig.  2.4  Upper  Section  of  the  Solid  Particle  Seeder. 
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Fig.  2.5  Lower  Section  of  the  Solid  Particle  Seeder. 

Worm  and  worm  gear  arrengement,  lead  screw, 
and  piston  are  shown. 


Fig.  2.6  Qualitative  Response  of  the  Seeder 
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is  distributed  evenly  and  injected  upward  inside  the  wind  tunnel  by  32 
equally  spaced  holes  of  about  4mm  in  diameter  just  upstream  of  the 
leading  edge  of  the  test  wall  boundary  layer.  The  basic  ideas  behind 
the  conical  manifold  are  to  construct  the  flow  paths  for  all  of  the  32 
injection  holes  equal  in  length,  and  to  eliminate  stagnation  points. 
This  provides  equal  tube  resistances  and  consequently  equal  flow  rates 
through  each  of  the  injection  parts. 

The  manifold  is  conical  in  shape  and  consists  of  thirty-two, 
1/4-inch  copper  tubes  inside  (Fig.  2.7).  The  spacings  between  the  tubes 
are  filled  with  sealing  foam  to  let  the  powder  flow  only  through  the 
tubes.  The  end  of  the  copper  tubes  are  connected  to  the  injection  ports 
via  plastic  tubes. 
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Fig.  2.7  Conical  Manifold(arrow)  for  Introducing  Seeded 
Air  into  Test  Section. 
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CHAPTER  III 
INSTRUMENTATION 


3.1  THEORY  OF  OPERATION: 

The  basic  principle  of  laser  doppler  velocimetry  is  shown  in 
Fig.  3.1.  The  laser  beam  is  split  into  the  two  equal  intensity  beams 
that  are  made  to  intersect  each  other  at  the  measurement  position.  The 
intersecting  beams  form  "stationary"  interference  fringes  inside  an 
ellipsoidal  measurement'  volume  at  the  focal  plane  of  the  lens.  To  have 
parallel  fringes,  the  intersection  point  should  occur  at  the  waist  of 
each  beam  since  this  is  the  location  where  plane  wave  fronts  exit  (Durst 
et  al. ,  1976). 

The  fringe  pattern  consists  of  dark  and  bright  regions  of  equal 
spacing.  The  spacing,  ax,  is  given  by: 


2  sin  (e/2) 


Where  e/2  is  the  half  angle  between  the  intersecting  beams  and  A 
is  the  wavelength  of  the  laser  light. 

If  a  particle  (1  to  4  micron  diameter)  passes  through  the 
measurement  volume,  as  it  is  shown  in  Fig.  3.1,  it  will  scatter  light 
which  can  be  detected  with  a  photo  detector.  The  Doppler  frequency 
(Durst  et  al , ,  1976)  of  the  scattered  light  is  proportional  to  the 
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Generated  by  a  Particle  Passing  through  a  Fringe 
Pattern  Formed  in  Measurement  Volume. 
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particle  velocity  which  is  given  by: 

f  =  sin  e/2  (3.2) 

X 

In  which  |U|  is  the  absolute  value  of  the  instantaneous  velocity  com¬ 
ponent  of  the  particle  in  the  direction  perpendicular  to  the  optical 
axis  {Fig.  3.1). 

There  are  various  time  and  frequency  domain  instruments  used  for 
measurement  of  Doppler  frequency  from  which  the  velocity  of  the  sus¬ 
pended  particles  in  one  particular  direction  can  be  calculated. 

Seeding  the  flow  with  the  so-called  "right"  material  is  one  of 
the  important  factors  in  laser  .  Doppler  anemometry.  From  a  fluid 
mechanic  point  of  view,  one  must  insure  that  the  particles  are  faith¬ 
fully  following  the  flow.  This  in  turn  means  that  the  density  of  the 
particles  should  be  as  close  as  possible  to  that  of  the  flowing  fluid. 
The  shape,  size,  and  concentration  of  the  particles  are  also  determining 
factors  for  their  motion. 

The  system  described  above  is  not  directionally  sensitive;  in 
other  words,  particles  moving  perpendicular  to  the  optical  axis  having 
the  same  speed  will  have  the  same  Doppler  frequency.  To  obtain  velocity 
directionality,  an  acousto-optic  modulator  called  the  Bragg  cell, 
(described  in  section  3.2.2)  is  used  to  modulate  the  incoming  laser 
beam. 

The  signal -to-noise  (SNR)  ratio  of  the  signal  output  of  the 
photo  detector  is  closely  related  to  the  ratio  of  the  particle  size  to 
fringe  spacing.  It  is  worth  mentioning  that  the  particle  size  normally 
refers  to  the  mean  particle  size  for  a  poly-dispersed  distribution.  The 
recommended  ratio  should  be  chosen  in  the  order  of  one  half.  Other 
factors  influencing  SNR  are:  Laser  power,  relative  intensities  of  the 
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two  light  teems,  alignment  of  the  tr*.  remitting  optics,  1  ight  collecting 
system,  and  clean  optics. 

Another  factor  worth  mentioning  is  particle  concentration.  If 
there  is  more  than  one  particle  at  c  time  inside  the  measurement  volume, 
constructi ve  or  destructive  superposition  of  the  scattered  lights  and 

consequently  of  electrical  detected  signals  can  occur  (Durst  et__al^ 

1976). 

3.2  _ Transmitt ing  Opt i cs 

Figure  3.2  shows  a  schematic  diagram  of  the  transmitting  optics 
looking  from  above.  An  Argon  Ion,  5-watt,  Spectra  Physics,  Model  104 
Laser  with  514.5  nm  wavelength  is  used  as  a  source  of  the  coherent, 
vertically  polarized  light  for  velocity  measurements.  The  Laser  beam  is 
passed  through  the  lens  combination  (LI  and  L2)  to  decrease  its  beam 
diameter  at  the  measruement  volume,  and  is  split  and  frequency  shifted 

by  a  Bragg  cell  to  make  the  system  directionally  sensitive.  The 

outcoming  beams  from  the  Bragg  cell  are  directed  towards  a  rotationally 
oscillating  scanner  mirror  (SM)  by  a  small  mirror  Ml.  The  beams  are 
then  reflected  to  the  last-stage  mirrors  (M5-M7)  by  a  long  and  fixed 
mirror  M3.  Measurement  volume  moves  perpendicular  to  the  plane  of  Fig. 
3.2  as  scanner  mirror,  SM,  oscillates. 

Requirements  for  the  components  in  the  transmitting  optics  are 
discussed  below. 

3 . 2.1 _ Lens  and  Mirrors 

A  finer  spatial  resolution  and  better  fringe  visibility  are 
usually  achieved  by  making  the  measurement  volume  as  small  as  possible. 
In  this  work,  there  is  a  minimum  limit  for  the  beam  diameter  at  the 
measurement  point.  This  limit  is  practically  determined  by  the  flatness 
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Table  3.1  Characteristic  Dimensions  of  the 
Transmitting  and  Receiving  Optics 


Components 

Parameter 

Value 

Mirrors  M3-M6 

length 

254  mm 

Lens  LI 

focal  length 

55  mm 

Lens  L2 

focal  length 
distance  between 

100  mm 

LI  and  L2 

160  mm 

Bragg  Cell  BC 

transducer  diameter 

25  mm 

transducer  frequency 

15.2  MHz 

Scanner  Mirror  SM 

size 

vertical  dimension  of 

26  mm  x  39  mm 

measurement  volumd 

^2  mm 

fringe  spacing 

*3.1  urn 

angle  between  intersecting 

beams 

distance  from  Laser  to 

the 

9.5° 

measurement  volume 

3708  mm 

Cylindrical  Lens  CL1 

focal  length 

300  mm 

size 

distance  from 

60  mm  x  200  mm 

measurement  volume 
to  CL1 

663  mm 

Cylindrical  Lens  CL2 

focal  length 

86  mm 

size 

57  mm  x  60  mm 

distance  between  Cl  &  i 

Cl 

546  mm 

distance  between  CL2  & 

PM  tube 

102  mm 

Mirror  M8 

size 

51  mm  x  381  mm 

Slit 

width 

1 . 5  mm 

length 

25.4  mm 
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of  the  mirrors  used. 

Two  circular  positive  lenses  are  used  after  the  laser  to 
decrease  the  beams'  diameters  at  the  measurement  position  down  to  2mm. 
The  estimate  of  the  beam  diameter  is  determined  by  intercepting  the  beam 
with  a  white  paper  at  the  measurement  point,  and  then  measuring  the 
diameter  of  the  circular  boundary  of  the  bright  spot  on  the  paper  as 
observed  by  a  special  laser  goggle. 

The  quality  of  the  mirrors  is  one  important  factor  in  this 
scanning  system.  For  instance,  scratches  on  the  mirror  could  increase 
the  noise  level  and  therefore  reduce  the  SNR.  A  long  and  flat  front 
surface  mirror,  which  is  aluminized  on  the  surface  nearest  the  object  to 
be  viewed,  is  used  so  that  light  reflects  from  this  near  surface  with  no 
refraction.  Mirrors  M3  -  M7  are  also  made  of  this  quality  (Fig.  3.2) 
and  are  mounted  on  aluminum  angles.  The  process  of  attaching  the 
mirrors  to  the  angles  and  the  technical  difficulties  observed  will  be 
discussed  later  in  this  section. 

The  last-stage  mirrors  consist  of  one  mount  with  two  mirrors 
(named  square  mirror),  M4  and  M5,  and  two  independent  mirrors,  M6  and 
M7.  They  are  all  mounted  on  a  separate  optical  table  which  is  fixed  to 
th'  traversing  table  (Fig.  3.4).  For  ease  of  alignment,  mirrors  M6  and 
M7  are  each  mounted  on  a  two-point  mirror  holder  that  can  rotate  the 
mirror  around  vertical  and  horizontal  axes.  The  square  mirror  has  three 
translational  and  one  rotational  degrees  of  freedom  which  make  the 
alignment  process  easy  to  accomplish. 

The  purpose  of  these  last-stage  mirrors  is  to  enable  one  to  vary 
the  fringe  spacing,  see  equation  2.1,  so  as  to  obtain  the  best  SNR  for 
the  seeding  used  here.  The  mirrors  M3  to  M7  should  be  long  enough  to 
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cover  the  desired  range  of  movements  of  the  measurement  volume. 

It  is  not  difficult  to  see  that  when  the  scanning  ranges  of  the 
two  beams  are  equal  in  length  on  the  paper  target,  located  exactly  at 
the  position  of  the  measurement  volume,  then  the  paths  travelled  by  the 
two  beams  are  equal  for  any  position  of  scanning  (path  length 
equalization  for  LDA  system).  In  fact,  this  is  one  of  the  advantages  of 
this  design  since  for  90%  of  the  cases  the  beams  can  be  aligned  to  scan 
the  same  length  by  only  adjusting  the  mirrors  M6  and  M7. 

The  main  difficulty  with  this  system  was  observed  when  the  two 
beams  (reflecting  from  mirrors  M6  and  M7)  were  aligned  to  cross  at  the 
measurement  point.  One  should  insure  that  these  two  laser  beams  inter¬ 
sect  with  each  other  and  are  perfectly  aligned  at  any  arbitrary  posi¬ 
tion  within  the  scan  range  (after  the  path  length  equalization  described 
above  has  been  accompl ished) .  For  a  simple  alignment  check,  the  system 
was  set  on  manual  scanning  mode,  and  it  was  observed  that  the  two 
"spots"  on  a  paper  target  (located  at  the  measurement  point)  were  not 
aligned  properly  as  shown  in  Fig.  3.5(b).  This  caused  significant  SNR 
reduction.  In  search  of  the  cause  of  this  problem,  flatness  of  the 
mirrors  was  suspected.  Therefore,  some  reliable  flatness  test  was 
needed.  For  a  simple  flatness  test,  a  horizontal  laser  beam  was 
reflected  from  various  points  on  each  mirror  to  a  graph  paper  target  a 
long  distance  away  as  the  mirror  mount  was  precisely  traversed  in  a 
vertical  direction  and  the  beam  displacement  was  observed  at  position 
(B)  as  it  is  shown  in  Fig.  3.6. 

Since  silicone  glue  was  originally  used  to  attach  the  mirrors  to 
the  angles,  this  technique  was  thought  to  disturb  the  beams'  alignment. 
Figure  3.7  shows  results  of  flatness  tests  conducted  on  two  mirrors 
before  and  after  gluing.  As  can  be  seen,  the  unit  vector  normal  to  the 
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Fig.  3.5  Alignement  of  two  Intersecting  Beams  as  Observed 

on  a  Paper  Target  Located  at  the  Measurement  Volume 
Position. 

(a)  Beams  are  separated. 

(b)  Beams  overlap. 

(c)  Perfect  alignement. 


Fig.  3.6  Optical  Arrangement  for  Mirror  Flatness 
(Surface  Curvature)  Test. 
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surface  of  the  mirror  M6  changes  drastically  as  one  scans  using  this 
mirror.  This  important  problem  was  solved  by  suspending  the  mirrors 
from  the  top  to  the  aluminum  angles,  using  a  simple  adhesive  tape,  which 
also  eliminates  any  changes  in  flatness  due  to  temperature  variation. 
The  results  of  flatness  test  conducted  on  the  existing  mirrors  are  shown 
in  Fig.  3.8. 

As  this  system  is  designed,  the  measurement  volume  scans  a  line 
rather  than  a  curve;  and  the  variation  of  the  fringe  spacing  is  at  most 
.8%  of  the  fringe  spacing  measured  when  the  two  intersecting  beams  are 
in  the  horizontal  plane,  which  was  found  to  be  immaterial  for  all 
practical  purposes. 

3.2.2  Bragg  Cell 

The  Bragg  cell  works  on  phenomena  known  as  Debye-Sears  effect 
and  Bragg  reflection.  Basically,  if  light  enters  a  material  with  an 
index  of  refraction  of  n,  there  will  be  a  phase  shift  occurring  on  the 
outcoming  light  from  the  material  which  depends  on  the  propagation 
constant  and  the  length  of  material.  If  this  n  can  be  changed  sinuso¬ 
idally  with  some  frequency,  f^,  then  the  outcoming  light  would  be  phase 
modulated  at  f  (Debye-Sears  effect). 

Debye-Sears  effect  in  Bragg  cell  is  observed  by  using  a  trans¬ 
ducer  vibrating  at  high  frequency  in  a  cavity  containing  distilled  water 
(this  creates  a  variable  index  of  refraction,  Fig.  3.9).  The  incoming 
laser  light  is  frequently  shifted  as  it  leaves  the  Bragg  cell.  The 
power  to  each  of  the  side  band  frequencies  can  be  changed  by  rotating 
the  Bragg  cell  around  the  axis  normal  to  the  paper  in  Fig.  3.9  (Bragg 
reflection). 


OBSERVED  DISPLACEMENT  WITH  RESPECT  TO 
ONE  END  OF  THE  MIRROR. 


ill 


5  10  15 

POSITION  01  W  ttiPJOR  ,  Ci. 


3.8  Chracteristics  of  the  Mirrors  Used  on 
the  Transmitting  Optics. 

Each  division  on  the  vertical  axis  is 
one  centimeter. 

▼  M6;0  M4;  ♦  M5;  A  M3;H  M7. 


Fig.  3.9  The  Debye-Sears  Effect. 
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The  Bragg  cell  transducer  used  in  the  present  system  has  a 
nominal  frequency  of  5.042  MHz  and  is  polished  for  the  third  overtone. 
It  is  important  to  drive  it  at  the  natural  frequency  so  that  maximum 
input  power  can  be  absorbed  and  the  gold-plated  surface  of  the 
transducer  will  not  peel  off.  The  natural  frequency  of  the  transducer 
was  measured  to  be  5.067  MHz. 

The  driver  unit  of  the  Bragg  cell  transducer  consists  of 
Heathkit  C.  B.  radio  transmitter.  Model  DX-60,  having  a  crystal 
oscillator  of  7.600  MHz.  It  amplifies  the  second  harmonic  to  give 
15.200  MHz  output.  The  output  of  the  transmitter  was  matched  to  the 
input  impedance  of  the  transducer  via  an  L-C  circuit  to  prevent  power 
reflection  back  to  the  transmitter. 

Using  Bragg  cell,  one  can  achieve  directional  sensitivity  for  a 
LDA  system.  In  the  present  system  the  first  order  frequency  shifted  and 
unshifted  beams  are  used  as  intersecting  beams.  These  two  beams  create 
a  fringe  pattern  that  looks  as  though  it  were  moving  with  velocity  equal 
to  Ax.  f  .  Here,  f  is  the  transducer  frequency,  and  Ax  is  fringe 
spacing. 

In  the  existing  system  fringes  are  moving  downstream  so  that  the 
modified  equation  for  the  Doppler  frequency  is: 

f  =  fm  -  (2U/A)  sin  8/2  (3.3) 

Where  f  is  Bragg  frequency,  U  is  the  instantaneous  velocity  normal  to 
the  optical  axis  of  the  system  (U  >  0  for  forward  flow  and  U  <  0  for 
blackflow).  The  Bragg  cell  is  covered  by  an  aluminum  housing  to  prevent 
rf  -  radiation  leakage  to  the  signal  processing  electronics. 
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3.2.3  Scanner: 

The  heart  of  the  scanning  LDA  is  a  6300PDT  scanner  from  General 
Scanning  Inc.  which  causes  the  measurement  volume  to  scan  through  the 
boundary  layer  in  the  test  section.  The  G300PDT  optical  scanner  is  a 
moving  iron  galvanometer  with  a  position  transducer  designed 
specifically  for  closed-loop  operation.  The  excursion  range  is  up  to  25 
degrees  of  mirror  rotation  (corresponding  to  50  degrees  of  optical  scan¬ 
ning).  The  position  transducer  operates  by  detection  of  capacitance 
variation  between  the  rotating  armature  and  a  set  of  stationary 
electrodes. 

Shaft  wobble  is  typically  below  five  arc-seconds.  The 
transducer  linearity  exceeds  ±  0.15%  of  the  working  excursion;  the  SNR 
affords  a  resolution  of  one  second  of  arc;  and  the  signal  response  time 
is  10  nsec. 

The  CCX102  controller  is  a  servo  control  driver  amplifiers  and 
includes  a  stepless  DC  heater  control  to  regulate  the  temperature  of 
G300PDT  scanner  (equipped  with  a  thermister  and  an  adhesive-bonded 
electrical  heater). 

Commercial  grade  scanner  mirror,  part  #  Ml-3226-39,  with 
thickness  of  3.2mm  and  dimension  of  26mm  x  39mm,  has  ALSiO  coating 
(aluminum  front  surface  over  coated  with  silicon  monoxide)  and  is  flat 
to  one  wavelength  per  centimeter. 

Since  the  signal  processing  described  in  Chapter  IV  depends  on 
the  linearity  of  the  ramp  section  of  the  scanner  internal  sawtooth 
generator  (which  drives  the  scanner  mirror),  a  test  is  conduced  to 
examine  its  behavior.  The  scanner  control  has  an  output  voltage  which 
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should  be  linearly  proportional  to  the  angular  position  of  its  mirror 
and  consequently  to  measurement  volume  position.  A  sewing  needle  is 
used  to  reflect  the  light  when  the  measurement  volume  passes  the 
location  of  this  needle,  and  a  photomultiplier  with  a  pinhole  detects 
this  light.  The  output  of  the  PM-tube  is  a  pulse  which  occurs  at  the 
time  when  the  measurement  volume  crosses  the  needle.  The  pulse  width 
depends  on  the  speed  of  the  scanning  as  well  as  the  needle  diameter. 
Fig.  3.10. 

The  position-voltage  output  of  the  scanner  control  unit  and  the 
output  of  the  PM-tube  are  fed  into  a  sample-and-hold  circuit  to  sample 
and  hold  the  position  voltage  when  the  measurement  volume  crosses  the 
needle  (see  Fig.  3.10).  The  position  of  the  needle  with  respect  to  some 
reference  level  is  accurately  determined  by  a  cathetometer. 

Figure  3.12  shows  the  scanner  position  voltage  output  signal 
versus  position  of  the  measurement  volume  (needle)  for  two  different 
frequencies.  The  correlation  coefficient  for  the  least  square  fit  was 
found  to  be  about  .999956,  which  shows  the  degree  of  linearity  of  the 
waveform.  A  convenient  calibration  technique  using  an  HP  9825B  computer 
was  used  in  practice  which  is  described  in  Chapter  IV. 

The  scanner  drift  was  determined  by  an  experiment  to  be 
approximately  0.0002  volts/min,  and  its  effects  on  the  calibration  curve 
are  practically  negligible. 
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Floor 


TOP  VIEW 


Fig.  3.10  Arrangement  for  Testing  the  Linearity 
of  the  Scanner. 


Fig.  3.11  Output  of  the  PM-Tube  and  Scanner  Position 

Voltage  Signal  as  Observed  on  a  Dual  Channel  _ Scope . 
Mote  the  voltage  sampled  by  Sample&Hold  in  Fig.  3.10. 


Sampled  Scanner  Position  Voltage  vs.  Position 
of  the  Needle(f1easurement  Volume). 
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3.3  Receivi  n  g_opjt  i  cs 

The  main  idea  in  the  design  of  the  receiving  optics  is  to  enable 
one  to  detect  the  scattered  signal  from  the  seed  particles  for  any 
position  of  the  measurement  volume  throughout  the  scan  range. 

Since  there  is  access  to  a  forward-scattering  signal  on  the 
opposite  side  of  the  wind  tunnel,  and  the  design  is  such  that  no 
simultaneous  movement  of  incident  and  receiving  optics  is  required, 
there  is  no  advantage  in  working  with  back  scattered  signal.  Figure 
3.13  shows  the  receiving  optics  arrangement.  The  assembly  is  isolated 
optically  inside  a  black  box  having  a  side  door  for  any  adjustment. 

3.3.1  Schematic  diagram  and  description 

Scattered  light  enters  into  the  box  through  four  cylindrical 
lens  (CL1),  50  mm  x  60  mm,  having  a  focal  length  of  300  mm.  The  choice 
of  this  focal  length  was  based  on  the  facts  that  the  distance  between 
the  center  line  of  the  tunnel  to  the  side  wall,  where  the  scattered 
light  is  gathered,  is  approximately  457  mm,  and  the  enlargement  of  the 
image  of  the  measurement  volume  in  the  streamwise  direction  on  the 
PM-tube  was  desired. 

To  make  the  length  of  the  receiving  optics  as  compact  as 
possible,  a  long  and  wide  front  surface  mirror,  M8,  of  the  same  quality 
as  mentioned  for  transmitting  optics  was  used  to  direct  the  scattered 
1  ight  downward. 

The  purpose  of  the  second  cylindrical  lense  (CL2),  57  mm  x  60  mm 
and  focal  length  of  86  mm,  is  to  construct  the  image  of  the  measurement 
volume  on  the  PM-tube  for  any  positions  of  scanning.  A  5  cm  long 
adjustable  slit  is  made  of  two  pieces  of  0.005  inch  thick  shim  stock  to 
gather  signal  from  a  portion  of  the  measurement  volume. 


Fig.  3.13  Cylindrical  Lens  Conf iguration  Schematic 
Diagram. 

See  Table  3.1. 
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Mirror  M8,  second  cylindrical  lens  (CL?),  and  PM-tube  are  all 
mounted  on  a  rack  such  that  they  can  be  moved  in  the  directions  shown  in 
Fig.  3.13.  There  are  two  long  slots  on  the  rack  so  that  the  holders 
can  be  fixed  at  any  vertical  position  desired.  To  have  additional 
flexibility  of  alignment,  the  rack  itself  can  be  moved  horizontally. 
The  first  cylindrical  lens,  CL1,  is  mounted  on  a  small  channel  and  its 
position  is  also  adjustable.  The  receiving  optics  arrangement  in  all  of 
the  experiments  reported  here  was  off-axis,  downstream  by  approximately 
25  degrees. 

Figure  3.14  shows  the  ray  tracing  for  the  receiving  cylindrical 
lenses.  Since  the  second  cylindrical  lens  only  works  in  a  vertical 
plane  one  can  use  a  simple  lens  equation. 

L_  +  L  =  ]_  (3.4) 

z+xx  X*  f2 

Where  Z  is  the  distance  between  two  cylindrical  lenses,  Xj  and  X£  are 
the  distances  from  the  measurement  volume  to  the  first  cylindrical  lens 
and  from  second  lens  to  the  PM-tube  respectively,  and  f^  is  the  focal 
lens  of  the  second  lens  (CL2). 

The  ratio  of  the  vertical  scan  range  to  the  image  length  on  the 

slit  is: 

Oj/If  =  (Z  +  Xx)/X£  (3.5) 

From  (3.4)  and  (3.5),  one  has: 

Xj  +  Z  =  1  +  (Oj/If )f2  (3.6) 

Xj  -  1  +  (lf/0j)f2  (3.7) 

ts&SiUSmc**  -■■■  ...  .  _  ....  - 
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Fig.  3.14  Ray  Paths  for  the  Receiving  Optics. 
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Fig.  3.15  The  Receiving  Optics  Assembly  inside  of  the 
Black  Box. 
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Using  the  same  lens  equation  for  the  first  cylindrical  lens 
which  works  in  horizontal  plane: 

1 


1  1 

—  +  — 


(3.8) 


X1  fl 


In  which  X  j  is  the  distance  from  CL1  to  the  PM-tube  and  f^  is  its  focal 
length.  From  (3.8)  finding  Xj  which  is  equal  to  Z  +  X£: 


Z  +  X0  = 


X1  '  fl 


2  Xj  -  fx 


(3.9) 


From  (3.6),  (3.7)  and  (3.9)  one  can  easily  derive  an  algebraic  equation 
which  determines  X^  for  any  choices  of  0^/ I ^  ratio  and  f^: 

f(x, )  =  Xn  +  (f1X1)/(X1-f1)  -  [2  +  (0,/If  +  VO,)]  f,  =  0  (3.10) 


Using  equation  3.10,  one  can  select  the  lenses  considering  the 
following  restrictions: 

1.  Since  the  photo  sensitive  front  surface  of  the  PM  tube 
is  almost  50  mm  in  diameter  and  the  image  of  the  scan 
range,  1^,  is  preferred  to  fall  on  the  central  region 
of  this  surface;  thus,  the  image  length  was  decided  to 
be  approximately  30  mm. 

2.  Since  the  maximum  boundary  layer  thickness  in  the  test 
section  is  roughly  330  mm,  the  maximum  scan  length, 

(0^),  was  chosen  to  be  381  mm. 

3.  The  distance  from  the  measurement  volume  to  the  first 


3. 
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cylindrical  length,  X^,  should  be  approximately 
greater  than  457mm  and  less  than  558mm,  since  the  half 
spanwise  width  of  the  tunnel  is  about  457  mm. 

Using  the  design  equation  (3.10)  with  f^  =  300  mm,  0^  =  381  mm 
and  1^  =  30  mm,  one  could  determine  the  focal  length  of  the  second 
cylindrical  lens  of  approximately  f„  =  86  mm  for  =  485  mm. 

Furthermore,  the  ratio  of  the  thickness  of  the  image,  d',  to  the 
streamwise  dimension  of  the  measurement  volume  would  be  about  1.62, 
allowing  easy  adjustment  of  the  slit  width  for  the  best  signals.  The 
width  of  the  first  cylindrical  lens  determines  the  solid  angle  from 
which  light  is  collected  and  should  be  matched  with  the  width  of  the 
second  lens. 

The  second  cylindrical  lens  should  be  wide  enough  to  gather  all 
the  scattered  light  that  is  passed  throughout  the  entire  length  of  the 
first  one.  The  existing  width  of  the  second  cylindrical  lens  is  57  mm, 
which  was  found  to  be  adequate  (see  Fig.  3.14). 

3.3. 2 _ Other  Alternatives  examined  for  Rece i ving  Optics. 

The  decision  on  the  design  explained  in  section  3.2  was  made 
based  on  examining  other  various  designs  of  which  two  are  worth  mention¬ 
ing: 

In  one  design  a  circular  lens  was  employed  to  construct  the 
image  of  the  probe  volume,  for  all  the  positions  of  scanning,  on  a  long 
vertical  slit  behind  which  the  central  region  of  a  fresnel  lens  was  used 
to  project  the  final  image  on  the  PM-tube.  It  was  found  that  this  final 
image  was  circular  in  shape  for  all  scan  positions.  It  is  worth  men¬ 
tioning  that  the  central  strip  of  the  fresnel  lens  works  basically  as  a 
very  wide  cylindrical  lens. 


ft*4 
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The  advantage  of  this  design  is  that  it  allows  one  to  use  a 
large  circular  lens  to  collect  more  scattered  light;  disadvantages  are 
the  poor  quality  of  the  fresnel  lens,  large  bulky  set  up,  and  poor 
signal-to-noise  ratio  at  the  ends  of  scanning. 

In  another  design  a  circular  lens  was  used  to  construct  the 
image  of  the  measurement  volume,  for  all  the  positions  of  scanning,  on  a 
short  slit  fixed  in  front  of  the  PM-tube  surface.  The  problem  of  poor 
SNR  on  the  scan  extremes  has  been  solved  by  this  design,  but  the  image 
on  the  PM-tube  was  very  narrow  and  consequently  very  sensitive  to  the 
slit  width.  The  only  advantages  of  this  design  are  the  compactness  of 
the  system  and  the  use  of  a  circular  lens. 

The  decision  on  the  final  design  used  in  this  work  was  made  on 
the  SNR  observed  by  spectrum  analyzer  (using  DOP),  and  the  fact  that  the 
counter  type  processor  needs  high  SNR  to  give  any  faithful  readings 
(Adrian  et  al .  1975). 

3.4  Fringe  Analysis  for  Scanning  LDA 

A  simple  equation  can  be  derived  to  represent  the  number  of 
fringes  crossed  by  a  particle  as  it  enters  the  measurement  volume.  The 
analysis  is  based  on  the  assumption  that  the  probe  volume  is 
two-dimensional  and  circular  in  shape,  and  the  particle  motion  is  domin¬ 
ated  by  its  streamwise  velocity  Up^.  In  a  directionally  sensitive  LDA 
system  the  fringe  pattern  is  moving  with  velocity  of: 

Uf  =  Ax  .  fm  (3.11) 

In  which  f  is  the  Bragg  cell  transducer  frequency  and  ax  is  fringe 
spacing.  In  the  system  employed  for  this  work,  fringes  move  in  the 
streamwise  direction  as  shown  in  Fig.  3.16.  Here,  the  relative  velocity 


Fig.  3. 


16  Number  of  Fringes  Crossed  by  a  Particle 
as  a  Function  of  IL /Vc  and  U  /Vc  'from 
Equ.  3.13.  f  S 

Arrows  on  each  curve  show  the  extremes 
of  the  U  /Vs  ratio  at  the  streamwise 
location  of  436.9  cm. 
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of  a  particle  with  respect  to  the  measurement  volume,  U  ,  is  referred  to 
when  the  probe  volume  scans  away  from  the  bottom  wall  of  the  wind  tunnel 
only.  Hence,  the  time  needed  for  a  particle  to  travel  one  fringe 
spacing  entering  from  the  position  angle  of  $  with  the  streamwise 
velocity  of  Up^  is: 

t  =  AX  (3.12) 


In  which  the  negative  value  for  Up^  represents  the  backflow  situation. 
The  resident  time  for  the  particle  in  the  probe  volume  is  (Fig.  3.16): 

y  AB  d  cos  (6  +  $-tt)  13j 

U  (U  2  +  V2)* 

r  '  px  s' 

Where  Vs  is  the  upward  scan  velocity  and  6  is  the  angle  between 
the  positive  direction  of  U  and  that  of  Ur  vector  (Fig.  3.16).  There¬ 
fore,  the  number  of  fringes  crossed  by  a  particle  is 


l  V  -  uf  i 


Ns  '  1  <v  +  vl> 


j~  (d/Ax)  cos  (e  +  <j>  -  n  ) 


(3.14) 


which  can  be  represented  in  non-dimensional  form  as: 


N$AX 

d 


i  i  -  <yv5xvs/upx) 

0  ♦  <w!]i 


cos(e  +  $  ) 


(3.15) 


i ,  i -v ■ 
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For  each  fixed  upswing  scan  velocity,  the  ratio  U^/V  is  a  constant 

number  and  the  right  hand  side  of  the  previous  equation  is  function  of 

V  /U  and  <t.  Note  that  the  Eq.  3.13  reduces  in  the  limit  to  the 
s  px 

stationary  measurement  volume  as  Vs/UpX  goes  to  zero. 

From  a  practical  point  of  view,  the  maximum  number  of  fringes 

crossed  by  a  particle  is  of  interest  and  can  be  derived  by  setting  the 

cosine  term  in  Eq.  3.15  equal  to  -1,  which  physically  means  that  the 

particle  passes  through  the  center  of  the  measurement  volume.  Figure 

3.16  shows  the  maximum  value  of  N  AX/d  versus  U  /V  for  three  values 

s  px  s 

of  Uf/V  ratio.  The  arrows  on  each  curve  show  the  extremes  of  the 

UpX/V$  parameter  at  the  streamwise  location  of  436.9  cm  in  the  wind 
tunnel.  The  maximum  number  of  fringes  crossed  by  a  particle  decreases 
as  the  scan  speed  is  increased  (low  values  of  U^/Vs). 

The  maximum  number  of  cycles  per  signal  burst  used  by  the 

counter  type  processor  in  this  work  is  32.  This  si.nplified  model  was 

used  to  insure  that  there  was  a  sufficient  number  of  fringes  crossed  by 
a  particle  for  various  scan  speeds  and  flow  velocities  in  the  boundary 
layer.  For  example,  with  a  maximum  streamwise  velocity  of  17  m/sec,  a 
fringe  spacing  of  3.1  pm,  a  fringe  velocity  of  47  m/sec,  and  a  small 
effective  streamwise  measurement  volume  dimension  of  3/4  mm  viewed  by 
the  slit;  a  particle  will  cross  no  more  than  32  fringes  with  an  upper 
limit  scan  speed  of  approximately  200  m/sec,  which  is  much  greater  than 
used  here.  To  use  a  higher  scan  velocity  and  still  obtain  a  sufficient 
number  of  fringe  crossings,  a  higher  Bragg  cell  frequency  could  be  used. 
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3.5  Velocity  Bias  Consideration 

According  to  velocity  biasing  analysis  of  McLaughlin  and 
Tiederman  (1973),  the  data  rate  in  a  "pointwise"  LDA  system  is 
proportional  to  the  magnitude  of  the  instantaneous  velocity  vector  U. 
Therefore,  the  velocity  amplitude  probability  distribution  is  biased 
towards  higher  velocities.  Likewise,  the  data  rate  in  "scanning"  LDA 
system,  in  which  measurement  volume  itself  moves,  is  proportional  to  the 
magnitude  of  the  instantaneous  "relative"  velocity  of  the  fluid  with 
respect  to  the  measurement  volume  U  . 

Based  on  the  refined  biasing  analysis  of  Buchhave  et  al .  (1979), 
the  percentage  biasing  error  decreases  as  the  parameter  Q  and/or  o/O 
decreases.  The  Q  parameter  is  the  ratio  of  the  minimum  number  of  zero 
crossings  required  for  a  valid  measurement  to  the  maximum  number  of 
fringe  crossings  available,  and  o/O  is  the  ratio  of  the  RMS  of  the 
velocity  fluctuations  to  the  mean  value. 

One  can  use  the  results  of  the  percentage  biasing  error 
mentioned  in  Buchhave  et  al . ( 1979)  to  have  an  estimate  of  the  biasing 
error  in  the  present  scanning  system,  but  one  should  note  that  the  mean 
"relative"  velocity  of  the  fluid  with  respect  to  the  measurement  volume, 
U^,  should  be  used  instead  of  the  mean  flow  velocity,  0,  in  the 
a/0  parameter.  For  a  given  o,  using  Fig.  3-5  of  Buchhave  et  _a_l_.  (1979), 
one  can  see  that  the  biasing  error  decreases  using  scanning  LDA,  in 
comparison  with  the  pointwise  LDA,  considering  the  fact  that  the  mean 
relative  velocity  of  the  fluid  with  respect  to  measurement  volume 
(scanning)  is  always  greater  than  the  mean  velocity  of  the  fluid 
(pointwise),  0.  It  is  interesting  to  note  that  the  biasing  error 
decreases  as  the  scan  velocity  is  increased. 
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Finally,  since  the  required  number  of  zero  crossings  for  a 
frequency  measurement  is  32,  and  the  effective  total  number  of  fringes 
available  is  about  300,  the  Q  factor  mentioned  above  is  almost  0.1,  and, 
therefore,  the  directional  dependence  is  very  weak  in  this  work. 
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CHAPTER  IV 

DATA  ACQUISITION  AND  PROCESSING 

Ideally,  in  a  scanning  laser  anemometer  system  one  needs  to  know 
the  exact  position  of  the  measurement  volume  at  the  time  when  a  "vali¬ 
dated"  velocity  measurement  is  detected.  The  data  acquisition  and 
processing  described  in  this  chapter  was  a  choice  based  on  available 
facilities  in  SMU.  It  might  not  be  the  best  approach  for  a  scanning  LDV 
system,  but  it  was  found  to  be  quite  adequate  (Durst  et  aK_  1981). 

4.1  General  description 

As  mentioned,  the  "ideal"  technique  is  to  sample  both  position 
voltage  output  of  the  scanner  driver  (a  sawtooth  waveform)  and  the  valid 
frequency  (or  velocity,  see  Eq.  3.3)  measurements  from  the  analog  output 
of  the  counter  processor  simultaneously.  The  whole  process  could  be 
triggered  via  "data  ready"  output  of  the  counter  which  sends  a  pulse 
every  time  a  new  measurement  has  been  validated. 

Since  the  HP-9825B  computer  used  in  this  work  cannot  be  pro¬ 
grammed  in  assembly  language  and  HPL  language  is  not  fast  enough  to 
switch  and  sample  position  as  well  as  counter  output  from  two  A/D 
convertors,  a  different  technique  for  data  acquisition  than  the  "ideal" 
one  was  used. 

It  is  worth  mentioning  that  the  switching  process  from  one  A/D 
converter  to  the  other  using  a  HPL  command  takes  on  the  order  of  5  msec. 
With  15  Hz  scanning,  using  scanner  internal  sawtooth  generator,  the  scan 
velocity  for  the  ramp  part  of  the  signal  is  about  5  m/sec.  Hence, 
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during  this  switching  time  the  measurement  volume  will  translate  by 
about  25  mm.  This  causes  two  problems:  firstly,  it  practically  elimi¬ 
nates  the  validated  signals  during  the  switching  time;  secondly,  posi¬ 
tion  ambiguity  occurs. 

The  retrace  part  of  the  sawtooth  position  waveform  was  not  used 
for  processing  due  to  its  poor  linearity.  The  symmetrical  waveform  is 
probably  preferred  for  the  external  control  of  the  scanner  driver,  since 
one  can  acquire  validated  signals  from  the  counter  processor,  both  when 
measurement  volume  scans  upward  and  downward. 

Considering  the  aforementioned  problems,  the  counter  processor 
analog  output  has  been  added  to  a  "marker"  signal  (stretched  SYNC  pulse 
in  Figs.  4.7  and  4.10),  which  marks  the  end  of  each  upswing  scan  and  was 
fed  into  a  fast  analog-to-digital  (A/D)  converter.  The  computer  can 
acquire  this  signal  (output  of  A/D)  as  fast  as  10,000  samples  a  second. 

It  should  be  emphasized  that  the  computer  does  not  sample  each 
time  a  validated  counter  measurement  is  available,  and  for  this  reason, 
some  memory  locations  will  be  wasted.  Data  processing  relies  heavily  on 
the  linearity  and  repeatability  of  the  upsweep  section  of  the  internal 
sawtooth  generator  used  by  the  scanner  drive  control  unit  (Brosens, 
1976).  This  will  be  discussed  in  more  detail  later  in  this  chapter. 

4. 2  Counter  Type  Processor 

As  mentioned  in  section  3.1,  the  scattered  light  from  particles 
passing  through  the  measurement  volume  is  detected  by  a  photomultiplier 
which  will  look  like  Fig.  4.1.  The  low  frequency  variation  corresponds 
to  the  passage  of  each  particle  through  one  or  other  or  both  intersect¬ 
ing  beams,  and  the  high  frequency  represents  the  passage  of  a  particle 
through  the  fringes. 
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The  TSI  Inc.  Counter,  Model  No.  1980,  is  used  for  finding  the 
high  frequency  component  of  the  signal  (Doppler  frequency)  which  is 
proportional  to  the  velocity  of  the  particle  (Eq.  3.1).  The  signal 
conditioning  section  of  the  TSI  counter  consists  of  a  signal  amplitude 
limit  control  which  eliminates  the  large  amplitude  signals  generated  by 
large  particles  which  are  not  following  the  flow,  low  and  high  pass 
filter  to  remove  the  low  frequency  variation  and  noise  respectively. 
Fig.  4.1,  and  finally  an  amplifier  with  adjustable  gain. 

Zero  crossing  technique  is  the  basic  principle  of  operation  used 
for  finding  the  Doppler  frequency.  The  zero  cross  voltage  level  is 
fixed  to  +50  mv  and  the  signal  amplitude  can  be  varied  by  the  gain 
control  of  the  amplifier.  The  analog  output  of  the  counter  is  a  voltage 
proportional  to  the  Doppler  frequency  which  in  turn  is  related  to  flow 
velocity  through  Eq.  3.3.  This  analog  output  is  used  for  all  data 
acquisitions  in  this  work.  The  counter  holds  the  output  voltage  of  the 
validated  measurement  until  the  next  one  occurs,  and  it  measures  the 
time  for  N  cycles  and  compares  them  with  that  of  N/2  cycles.  A  division 
is  actually  performed  so  that  it  is  a  true  percentage  comparison  which 
can  be  set  via  a  dial  on  the  counter. 

Figures  4.2  (a)  and  4.2  (b)  represent  the  typical  signals 
observed  on  the  oscilloscope  using  baby  powder  as  seeding  agent  while 


scanning. 
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Fig.  4.1  Characteristics  of  Electrical  Signal. 

(a)  Signal  from  photomultiplier  tube. 

(b)  After  high-pass  filtering  to 
eliminate  dc  pedestal. 
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4 . 3_  Po  s  it  ion  Cali  bra  tion  of  the  Measurement  Volume 

One  needs  to  know  the  position  at  which  a  validated  measurement 
occurs.  For  this  purpose,  the  scanner  driver  has  a  position  voltage 
output  (sawtooth)  which  is  proportional  to  the  position  of  the  measure¬ 
ment  volume  and  can  be  used  for  calibration. 

A  special  black  paper  target  having  15  holes  of  1mm  diameter 
and  25.4  mm  apart  vertically  is  attached  to  a  stable  frame  such  that  the 
first  hole  is  25.4  mm  from  the  tunnel  floor.  A  white  piece  of  paper  is 
attached  behind  each  of  the  holes  to  scatter  light  as  laser  beams  pass 
over. 

The  calibration  procedure  starts  by  locating  the  paper  target  in 
the  position  where  the  velocity  profile  is  needed  such  that  the  measure¬ 
ment  volume  scans  through  the  holes  on  the  target.  The  same  receiving 
optic  arrangement  is  used  for  calibration,  and  as  the  measurement  volume 
passes  over  each  of  the  holes  a  pulse  is  generated  by  the  PM-tube. 
Since  the  pulses  on  the  retrace  are  not  desirable,  a  low  pass  filter  is 
used  to  attenuate  and  finally  eliminate  them  via  a  comparator  (Fig. 
4.3).  Hence,  only  pulses  corresponding  to  the  ramp  section  of  the  posi¬ 
tion  voltage  output  of  the  scanner  remain  after  the  comparator.  The 
combination  of  an  inverter,  a  comparator,  and  one-shot,  (see  section 
4.4.1)  make  the  pulses  look  like  Fig.  4.3  (G).  This  signal  is  used  to 
trigger  an  A/D  converter  externally  in  order  to  sample  the  position 
voltage  of  the  scanner.  The  HP  9825  computer  acquires  the  sampled 
position  voltages. 

The  software  written  for  the  calibration  acquires  50  cycles  of 
consecutive  scans  of  the  measurement  volume  and  computes  the  linear 
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Fig.  4.4  Effect  of  the  Scan  Range  Shift(c)  on  the 

Waveform(a)  and  on  the  Calibration  Line(b). 
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least  square  fit  through  the  sampled  position  voltages.  Therefore,  one 
has  a  linear  relation  between  the  position  voltage  output  of  the  scanner 
and  the  position  of  the  measurement  volume. 

It  was  found  that  the  scan  range,  once  set  by  the  scanner 
driver,  remained  constant;  but  the  whole  scan  range  shifted  either  up  or 
down.  Fig.  4.4  (c).  As  far  as  calibration  is  concerned,  this  scanner 
"shift"  does  not  have  significant  effects  either  on  slope  or  on  inter¬ 
cept  of  the  calibration  line  and  has  been  confirmed  experimentally.  It 
will,  however,  cause  problems  in  determining  the  position  corresponding 
to  each  computer-sampled  point  which  is  discussed  in  section  4.4.2. 

4.4  Data  Acquisition 

This  section  is  divided  into  two  parts.  The  first  part  is 
devoted  to  the  description  of  the  special  circuitry  made  for  the  acqui¬ 
sition,  called  "pulse  marker",  and  the  second  part  deals  with  the  data 
acquisition  logic  and  procedure. 

For  more  clarity  in  the  logic  and  procedure  of  the  acquisition 
and  processing,  the  reader  is  encouraged  to  refer  to  the  key  Figs.  4.7 
and  4.10  frequently. 

4.4.1  Pulse  Marker 

The  circuit  diagram  of  the  pulse  marker  is  shown  in  Fig.  4.5. 
A  dual  amplifier  is  used  as  voltage  follower  to  have  high  input  imped¬ 
ance  and  to  prevent  loading  of  the  equipment  from  which  signals  are 
taken. 

As  can  be  seen,  the  diagram  has  basically  two  branches  which  are 
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added  together  to  construct  the  output.  Typical  operation  of  the  pulse 
marker  which  is  closely  related  to  its  actual  one  in  data  acquisition  is 
explained  in  Fig.  4.6. 

Signal  "A"  in  Fig.  4.6  is  available  by  the  scanner  drive  control 
called  SYNC  pulse  and  always  occurs  at  the  sane  location  shown  in  that 
figure.  The  purpose  of  the  lower  branch  in  the  pulse  marker  circuit 
diagram  is  to  generate  a  variable  pulse  width  controlled  via  a  capacitor 
box  (Fig.  4.6  b). 

In  data  acquisition,  the  input  signal  (c)  in  Fig.  4.6  is  the 
analog  output  of  the  counter  processor.  The  voltage  variation  of  the 
counter  analog  output  for  the  settings  used  in  the  data  acquisition  was 
within  3.5  ±  1  volts.  Since  better  voltage  resolution  (.0005  volts)  was 
available  working  on  ±  1  volt  operation  range  of  the  fast  A/D  converter, 
the  upper  branch  of  the  pulse  marker  circuit  was  added  to  shift  the 
input  signal  (c)  by  an  appropriate  value,  E,  in  order  to  bound  this 
input  signal  within  ±1  volt  range. 

If  the  voltage  level  at  "A"  in  Fig.  4.5  exceeds  the  adjustable 
discrimination  voltage,  E^,  then  the  output  of  the  comparator  goes  high 
and  remains  there  and  firally  drops  to  zero  (low)  when  the  voltage  level 
at  "A"  goes  lower  than  Ed.  Theoretically,  as  far  as  data  acquisition  is 
concerned,  one  may  not  need  a  comparator;  but,  it  was  added  since  this 
same  circuit  was  used  for  calibration,  and  this  additional  feature  was 
found  to  be  practicable. 
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4.4.2  Data  Acquisition  Logic  and  Procedure 

The  whole  idea  is  that  if  one  knows  the  position  calibration 
line  and  the  sawtooth  position  waveform,  one  can  find  the  position 
corresponding  to  each  valid  velocity  measurement  in  this  scanning 
system.  Figures  4.7  and  4.10  explain  the  idea.  Step  by  step  details  of 
the  procedure  will  follow. 

As  mentioned  before,  data  acquisition  relies  heavily  upon  the 
linearity  and  repeatability  of  the  sawtooth  position  waveform  of  the 
scanner.  The  knowledge  of  this  waveform  is  needed  for  finding  the 
position  corresponding  to  each  piece  of  computer-stored  points  (output 
of  A/Dl  in  Fig.  4.7).  In  fact,  only  the  validated  velocity  measurements 
that  occur  during  upward  scan  of  the  measurement  volume  are  of  interest. 
Therefore,  software  was  written  to  sample  the  sawtooth  waveform  and  to 
compute  the  slope,  intercept,  minimum,  and  maximum  voltage  for  the 
upsweep  section  of  this  waveform  (statistical  analysis).  The  results  of 
this  analysis  showed  that  the  maximum  rms  of  the  voltage  variations  at  a 
given  phase  on  the  upsweep  part  did  not  exceed  0.2%  of  the  mean  voltage 
value  in  that  same  phase  which  represents  the  repeatability  of  the 
waveform. 

Although  the  shape  and  the  frequency  of  the  sawtooth  waveform 
are  adequately  stable,  but  the  whole  waveform  might  shift  up  or  down 
either  manually  or  by  itself  (Fig.  4.8),  what  it  basically  means  is  that 
the  intercept  for  the  upsweep  section  of  the  sawtooth  waveform  changes 
in  time.  The  knowledge  of  the  intercept  at  the  time  when  each  record  of 
the  frequency  data  is  taken  would  be  essentially  important  for  finding 
the  position  corresponding  to  each  frequency  data  point,  see  Figs.  4.8 


and  4.10. 
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The  SYNC  pulse  of  the  scanner  driver  can  be  used  to  keep  track 
of  the  amount  of  the  voltage  shift.  The  voltage  corresponding  to  the 
time  when  SYNC  pulse  occurs  is  called  "Reference  Voltage",  Fig.  4.8. 
From  the  statistical  analysis  of  the  sawtooth  waveform,  one  has  slope, 
intercept,  and  the  reference  voltage  at  the  time  this  waveform  was 
sampled  by  the  computer.  The  difference  between  the  reference  voltages 
for  each  record  of  the  frequency  data  acquisition  (shown  as  solid  line 
in  Fig.  4.8)  and  that  of  the  previously  known  value  when  the  sawtooth 
waveform  was  sampled  for  statistical  analysis  (broken  line.  Fig.  4.8) 
determines  the  correction  to  the  intercept  for  the  upsweep  part  of  the 
waveform. 

To  summarize,  the  first  step  in  the  data  acquisition  is  position 
versus  voltage  calibration;  the  second  step  is  the  statistical  analysis 
on  the  scanner  sawtooth  waveform;  and,  finally,  the  third  step  is  the 
acquisition  of  the  frequency  voltage  from  the  counter.  Figure  4.7  shows 
the  block  diagram  for  the  final  step  of  the  data  acquisition. 

In  Fig.  4.7  the  voltage  output  of  the  LDA  processor  (TSI  coun¬ 
ter),  which  is  proportional  to  the  Doppler  frequency  or  velocity  of  the 
particles  (fluid),  is  shifted  by  the  pulse  marker  (to  be  within  ±  1 
volt)  and  is  added  to  the  marker-pulse  (Figs.  4.6  and  4.10).  It  is 
important  that  the  pulse  width  of  the  marker-pulse  be  greater  than  the 
period  of  the  pulses  used  for  sampling  in  A/D  1  (section  4.5). 

The  final  phase  of  data  acquisition  starts  by  acquiring  100 
samples  from  the  reference  voltage,  just  before  the  frequency  data  is 
acquired,'  using  A/D2  with  SYNC  as  its  sampling  pulse  (Fig.  4.7),  then 
switching  to  the  A/Dl  to  acquire  a  single  record  of  17500  points,  with  a 
controllable  sampling  period  (pulse  generator)  of  T 


and,  finally. 


Pulse  marker  circuit  stretches  scanner  SYNC  pulse 
and  adds  analog  output  from  counter. 


Pig.  4.8  Shift  of  the  Sawtooth  Waveform. 
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switching  back  to  A/D2  and  to  acquire  another  sample  of  the  reference 
voltage,  see  Fig.  4.7. 

The  average  of  the  reference  voltages  sampled  in  the  beginning 
and  at  the  end  of  each  record  was  used  for  data  processing  and  monitor¬ 
ing  purposes.  The  difference  between  this  average  reference  voltage 
(when  this  record  of  frequency  data  is  taken)  and  that  known  when 
statistical  analysis  on  sawtooth  position  waveform  is  done  (in  second 
step  of  acquisition),  determines  the  intercept  correction  in  Figs.  4.8 
and  4.10.  Every  pulse  of  the  pulse  generator.  Fig.  4.7,  triggers  a  new 
conversion  cycle  for  the  fast  A/D  converter  which  lasts  50  psec.  Hence, 
the  period  of  the  pulses,  Tex^,  cannot  be  less  than  50  psec.  The  choice 
of  sampling  period  is  essential  both  in  the  acquisition  and  processing 
of  the  data.  Here,  this  choice  of  sampling  period  is  discussed  briefly 
in  regard  to  the  position  uncertainty  in  this  work. 

Figure  4.9  (a)  shows  a  short  portion  of  the  output  of  the  pulse 
marker  (see  Fig.  4.10  c)  which  has  been  digitized  by  the  A/Dl.  Pro¬ 
cessing  logic  in  this  work  accepts  the  sample  number  2  as  the  time  when 
particle  passes  the  measurement  volume,  while  in  reality  the  particle 
passed  the  volume  sometime  earlier.  The  maximum  error  in  time  when  the 
particle  crosses  the  intersection  point  of  the  two  laser  beams  is  equal 


If  the  measurement  volume  scans  a  length  equal  to  R  mm  with 
frequency  of  f  using  sawtooth  internal  generator  of  the  scanner,  then 
the  distance  travelled  by  the  measurement  volume,  aY^,  in  the  time 
interval  of  Tgxt  is  (Fig.  4.9  b/: 


sYr  ■  CR/(1-TR/TS)]  fs  Text  <») 


(4.1) 
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For  any  arbitrary  scan  range,  L-cm,  one  can  use  the  above  table  with  the 
correction  equation  AY.  =  aYr(L/R). 
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The  value  of  1  -  TR/TS  is  .921  for  15  Hz  frequency  of  scanning 
and  remains  fairly  constant  for  the  wide  range  of  frequencies  (Fig. 
4.9  b). 

Table  4.1  shows  the  distance  travelled  by  the  measurement  volume 
during  sampling  period  for  various  scan  frequencies  and  sampling  per¬ 
iods,  Tgxt  (R  is  taken  as  30.5  cm.).  Since  the  beam  diameter  is 
approximately  2  mm,  the  T^  for  each  scan  frequency  is  chosen  so  that 
the  position  uncertainty,  due  to  measurement  volume  movement,  becomes 
less  than  the  beam  diameter.  Hence,  before  the  beginning  of  the  final 
stage  of  data  acquisition,  one  should  make  an  appropriate  decision  of  a 
value  for  the  sampling  period.  There  is  another  important  determining 

factor  for  the  value  of  T  ,  which  will  be  discussed  in  section  4.5. 

ext 

4.5  Data  Processing 

Figure  4.10  shows  the  details  of  the  processing.  The  reference 
voltage  for  each  data  record,  section  4.4.2,  is  used  to  find  the  cor¬ 
rection  on  the  intercept  for  the  ramp  section  of  the  sawtooth  waveform 
for  each  particular  record  (Fig.  4.10).  Note  that  the  reference  of 
time  for  calculating  the  slope  and  intercept  for  the  upsweep  part  of  the 
waveform  is  taken  at  the  lowest  voltage  of  the  scan  on  the  ramp. 

The  height  of  the  marker-pulse  (stretched  SYNC  pulse.  Fig. 
4.10)  is  such  that  it  always  exceeds  the  saturation  voltage  level  (+1 
volts)  of  the  A/Dl.  This  pulse  marks  the  end  of  each  upswing  scan  cycle 
which  is  digitized  and  stored  by  the  computer.  For  this  reason,  the 


width  of  the  marker-pulse  must  be  slightly  greater  than  Text  to  insure 
that  the  A/Dl  samples  at  least  once  on  the  marker-pulse. 


Data  processing  program  starts  by  scanning  backward  in  time 
through  the  stored  record  searching  for  the  sampled  marker-pulses  and 
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Fig.  4.10  Data  Acquisition  and  Processing  Logic  and  Signals. 

(a)  Internal  sawtooth  waveform  of  scanner. 

(b)  SYNC  pulse  from  scanner. 

(c)  Sum  of  the  analog  output  from  counter  processor 
and  the  stretched  SYNC  pul se(Marker-Pul se) . 

(d)  Scanner  voltage  vs.  posi tion(Cal ibration  Line). 

(e)  Sampled  data  points  from(c)  with  corresponding 
position  voltage  on  the  linear  least-squre  fit 
waveform. 
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creating  a  new  record  by  eliminating  the  sample  points  on  the  retrace  of 
each  and  every  cycle  in  the  record.  In  this  new  record,  there  is  an 
indicator  for  the  beginning  of  each  upswing  scan  which  at  the  same  time 
represents  the  end  point  of  the  previous  scan  cycle,  see  Fig.  4.11. 

One  knows  the  corrected  intercept,  slope,  the  time  corresponding 
to  the  maximum  location  of  scanning  (t^  in  Fig.  4.10),  and  the  sampling 
period,  (T  t)  from  data  acquisition  process.  Therefore,  one  can 
determine  the  position  corresponding  to  each  computer-stored  point 
(output  of  A/D  1)  by  using  calibration  line  as  shown  in  Fig.  4.10. 

The  total  scan  length  is  divided  to  a  certain  number  of  position 
"bins".  The  processing  program  searches  through  the  computer-stored 
points  for  validated  velocity  measurements  in  each  bin  by  finding  the 
voltage  difference  of  the  two  consecutive  samples,  AE.  If  the  absolute 
value  of  this  difference  is  greater  than  a  certain  number,  then  the 
sample  which  occurred  later  in  time  is  taken  as  the  valid  velocity 
measurement  in  that  particular  bin.  Note  that  since  velocity  and 
Doppler  frequency  for  a  particle  are  related,  thus  the  terms  velocity 
measurement  or  frequency  measurement  are  equivalent  in  this  context. 

Since  the  voltage  resolution  of  the  A/Dl  is  .0005  volt,  the 
absolute  value  of  the  voltage  difference  (|aE|)  mentioned  above  should 
be  greater  than  this  number.  In  fact,  it  was  found  experimentally  that 
the  | AE |  must  be  greater  than  or  equal  to  0.005  volt  and  less  than  0.01, 
in  order  to  pick  the  true  validated  velocity  measurements.  Note  that 
the  0.005  voltage  difference  corresponds  to  approximately  ±  0.061  m/s, 
using  a  3.1  vm  fringe  spacing. 

At  the  final  stage  of  the  data  processing  program,  all  of  the 
validated  velocity  measurements  for  each  "position  bin"  are  used  in 
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computing  mean  and  RMS  values  for  that  particular  bin.  Furthermore, 
program  prints  out  the  number  of  valid  measurements  for  each  bin  which 
is  used  for  computing  the  data  rate  distribution  (chapter  V). 

From  the  data  processing  logic,  it  is  obvious  that  at  least  two 
computer-stored  points  are  needed  to  determine  a  single  validated 
measurement  in  each  position  bin  (see  Figs.  4.9  a  and  4.10  d).  If  one 
wants  to  have  greater  than  or  equal  to  "S"  number  of  computer  data 
samples  in  each  "position  bin,"  then  the  sampling  pulse  period,  T  t, 
must  be 


Text  <  (-«21/S){Ts/Mb> 


(4.2) 


where  Ts  is  the  period  of  sawtooth  scanning  waveform,  and  is  the 
total  number  of  position  bins  in  the  scan  range  (section  4.4.2).  So  the 
Text  should  be  chosen  such  that  to  minimize  the  position  uncertainly  as 
well  as  satisfy  equation  4.2. 

Table  4.2  represents  the  maximum  Text  allowable  in  order  to  have 
at  least  three  numbers  of  computer  data  samples  in  each  position  bin. 


4.6  Static  and  Dynamic  Check  of  the  Entire  System 

Due  to  the  lack  of  availability  of  the  appropriate  equipment, 
the  data  processing  as  explained  in  previous  sections  is  somewhat 
complicated.  Therefore,  a  test  procedure  was  needed  to  check  the  proper 
operation  of  the  whole  system  including  the  softwares. 

"*  There  are  two  tests  called  static  and  dynamic  checks.  The 
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Table  4.2  Maximum  sampling  pulse  period  allowable  in  order  to  have 
3  (=S)  computer  samples  in  each  position  bin  using  30 
(=Nj5)  number  of  position  bins  in  scan  range. 

VHz  Text->,sec 


6 

1700 

10 

1000 

15 

682 

20 

512 

25 

409 

30 

292 

40 

256 

50 

205 

60 

171 

70 

146 

100 

102 

f 
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static  test  has  a  two-fold  purpose  of  voltage  calibration  as  well  as 
testing  the  system.  The  terminology  of  "static"  and  "dynamic"  will  be 
explained  by  the  procedure  of  each  test. 

The  static  test  procedure  is  to  follow  all  the  steps  of  data 
acquisition  and  processing  but  using  a  constant  voltage  from  a  DC  power 
supply,  instead  of  the  TSI  counter,  as  shown  in  Fig.  4.12.  The  input 
voltage  is  measured  by  a  digital  voltmeter  (DVM). 

The  end  result  of  the  processing  is  the  mean  voltage  for  each 
position  bin  of  the  scanning  which  must  remain  the  same  for  all  position 
bins.  The  entire  procedure  could  be  repeated  for  various  fixed  input 
voltages  and  a  voltage  calibration  line  be  derived.  Figure  4.13  shows 
such  a  calibration  line  which  later  is  incorporated  in  the  data  process¬ 
ing  program. 

The  second  test  on  the  entire  system  uses  the  position  waveform 
of  the  scanner,  as  the  simulated  TSI  counter  output,  as  shown  in  F^'g. 
4.14;  and  for  that  reason  the  test  is  called  a  dynamic  test. 

As  the  data  processing  logic  works,  see  Fig.  4.10,  if  the 
software  is  to  follow  the  logic,  then  one  should  not  only  obtain  a  line, 
if  a  plot  of  the  mean  value  for  each  position  bin  versus  position  is 
constructed,  but  this  line  should  represent  the  ramp  section  of  the 
sawtooth  waveform. 

Figure  4.15  shows  the  results  of  such  a  test.  The  solid  line 
represents  the  result  taken  from  the  statistical  analysis  on  the  wave¬ 
form,  see  section  4.4.2. 

4.7  Position  and  Velocity  Uncertai nty 

In  sections  4.4.2  and  4.5  concerning  data  acquisition  logic  and 
processing,  it  was  referred  that  the  major  position  uncertainty  was  due 
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to  the  technique  used  for  determining  the  position  corresponding  to  each 
counter  validated  velocity  measurement  (Fig.  4.10). 

Table  4.1  shows  the  maximum  position  uncertainty  associated  with 
aforementioned  processing  technique  in  millimeter.  Since  the  processing 
program  accepts  the  point  No.  2  in  Fig.  4.9  (a)  as  the  time  when 
particle  passes  the  measurement  volume,  therefore,  one  has  the  true 
position  where  the  particle  crossed  the  fringes  as: 

(a,R 

'true  ■  V/2-  (4.3) 


Where  is  the  position  corresponding  to  the  point  No.  2  in  Fig. 
4.9(a),  d  is  the  beam  diameter,  and  aYr  is  given  in  Table  4.1.  Accord¬ 
ing  to  the  uncertainty  analysis  of  Kline  and  McClintock  (1953),  in  all 
of  the  experiments  in  this  work  the  Tgxt  in  Table  4.1  is  chosen  such 
that  the  maximum  position  uncertainty  due  to  measurement  volume  movement 
is  less  than  the  beam  diameter.  Hence,  the  dominant  uncertainty  would 
be  due  to  the  large  beam  diameter  itself. 

Another  source  of  position  uncertainty  is  due  to  the  fact  that  a 
linear  least  square  fit  to  the  ramp  upsweep  of  the  sawtooth  waveform  is 
used  which  even  for  high  scan  frequencies  (60Hz)  does  not  exceed  1  mm, 
and  therefore  is  negligible. 

The  uncertainty  due  to  velocity  is  mainly  associated  with  the 
counter  processor.  The  counter  measures  the  time  for  32  cycles  and 
compares  this  time  to  that  of  16  cycles.  The  comparison  is  done  on  the 
frequencies  calculated  based  on  32  and  16  cycles  by  actually  dividing 
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the  two  numbers.  The  comparison  limit  can  be  set  on  the  counter  from  1* 
to  20%  discrepancy  between  the  two  values.  The  equation  that  determines 
the  velocity  can  be  derived  using  equations  (2.1)  and  (3.1) 

(15.2  x  106  -f  •  106)  *  (Ax  •  10‘6)  =  U  (4.4) 

in  which  f  is  the  Doppler  frequency  measured  by  counter  in  MHz,  ax  is 
the  fringe  spacing  in  micron,  and  U  is  the  instantaneous  streamwise 
velocity  in  m/sec. 

In  equation  4.4,  an  error  in  the  measurement  of  the  Doppler 
frequency  of  ±df  is  equivalent  to  the  velocity  uncertainty  of  ±10.2  df 
using  fringe  spacing  of  3.1  ym.  For  1%  comparison  limit  the  ±df  is 
about  .150  MHz,  meaning  that  the  associated  velocity  uncertainty  is  at 
most  +.45  m/sec. 

The  TSI  counter  uses  a  250  MHz  clock  to  measure  the  frequency. 
Therefore,  uncertainty  due  to  +1  count  ambiguity.  Durst  et  at.  (1976), 
is  almost  ±.05  m/sec  for  maximum  forward  flow,  which  is  negligible 
compared  with  the  comparison  uncertainty  of  ±  .45  m/sec. 
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CHAPTER  V 

EXPERIMENTAL  RESULTS  AND  DISCUSSION 


5.1  Mean  and  RMS  Values  Using  Scanning  LDV 

There  were  several  questions  raised  concerning  the  operation  of 
this  new  design:  1.  Is  this  system  reliable  for  any  velocity  measure¬ 
ment?  2.  What  is  the  maximum  scan  frequency  achievable  in  this  design? 
3.  What  are  the  limitations?  Therefore  a  series  of  experiments  were 
made  in  an  attempt  to  answer  the  above  fundamental  questions. 

For  comparison  purposes,  in  all  of  the  experiments  the  results 
of  the  pointwise  measurements  using  the  same  scanning  LDA  setup  are 
compared  with  that  of  scanning  data  for  various  frequencies  in  a  turbu¬ 
lent  separated  flow  generated  in  the  SMU  wind  tunnel. 

Figures  5.1  and  5.2  show  mean  and  rms  velocity  profiles  for  each 
position  bin  and  for  various  scan  frequencies.  In  each  case  the  total 
record  consists  of  several  number  of  smaller  data  records  taken  during 
experiment;  the  total  record  time  in  previous  figures  for  each  frequency 
was  at  most  one  minute.  The  boundary  layer  thickness  for  Figs.  5.1  and 

5.2  is  about  35  cm.  Consequently  due  to  variation  of  the  vertical 
distance  of  the  turbulent-nonturbulent  interface  and  the  fact  that  much 
of  the  baby  powder  remains  inside  the  turbulent  boundary  layer,  the  mean 
values  are  biased  near  Y  s  25  cm,  besides,  the  particle  concentration 
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According  to  Yanta  and  Smith  (1973)  taking  a  flow  with  5%  level 
of  turbulence,  a  sample  size  of  166  is  necessary  to  obtain  a  .99% 
confidence  level  that  the  experimental  mean  differs  from  the  true 
population  mean  by  less  than  1%.  But  for  population's  standard  devia¬ 
tion,  3150  sample  size  is  required  for  a  confidence  level  of  95%  and  an 
error  of  2.5%.  The  fact  that  one  needs  larger  sample  size  for  standard 
deviation  than  the  mean,  can  also  be  seen  in  a  work  by  Bates  and  Hughes 
(1976),  in  which  they  concluded  that  for  sample  size  of  3500  and  above, 
the  turbulence  intensity  shows  no  variation  with  sample  size.  One  can 
see,  in  their  work,  that,  with  sample  size  of  500,  the  errors  for  mean 
and  turbulence  intensity  are  in  the  order  of  4%  and  22%  respectively. 

The  maximum  number  of  validated  counter  measurement  for  each 
position  bin  in  Figs.  5.1  and  5.2  does  not  exceed  350.  Near  the  wall 
and  edge  of  the  boundary  layer,  data  rates  are  low;  therefore,  one  has 
smaller  sample  size  for  corresponding  position  bins,  which,  in  turn, 
means  higher  percentage  error.  Despite  that,  good  agreement  exists 
between  the  mean  values.  Fig.  5.1,  up  to  a  scan  frequency  of  59  Hz  (scan 
speed  of  19  m/sec. ) 

One  should  consider  that  as  the  measurement  volume  scans  with 
this  rapid  velocity  through  the  highly  turbulent  separated  flow  the 
counter  processor  receives  variable  signal-to-noise  (SNR)  ratio  signals 
due  to  measurement  volume  movement  and  from  variations  in  seeding 
density  and/or  particle  size.  Orloff  et  al .  (1974)  also  reported  these 
effects. 

The  rms  of  velocity  fluctuations  shows  more  scatter,  but  is 
within  measurement  uncertainty  (section  4.7);  and  as  mentioned  before, 
one  reason  is  due  to  insufficient  number  of  validated  counter  measure¬ 
ments  in  each  bin  (small  sample  size).  This  is  especially  seen  for  the 
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position  bins  corresponding  to  the  extremes  of  the  scanning  range. 
Better  results  were  taken  by  proper  focusing  of  the  image  of  the  mea¬ 
surement  volume  on  PM-tube  and  increasing  the  sample  size  (Figs.  5.3  and 
5.4). 

Once  the  proper  operation  of  this  design  was  established,  the 
system  was  used  to  survey  the  flow  in  the  SMU  wind  tunnel.  Figures  5.3 
and  5.4  show  the  mean  and  rms  of  the  velocity  fluctuations  at  four 
different  streamwise  locations  in  the  wind  tunnel,  using  scanning  LDA, 
as  compared  with  those  of  pointwise  velocity  measurements  (using  the 
same  scanning  arrangement).  The  sample  size  for  each  position  bin  does 
not  exceed  600,  and  for  this  reason,  all  of  the  pointwise  measurements 
are  with  sample  size  of  500  to  compare  with  those  of  scanning  data. 

Figure  5.5  shows  the  intermittency  factor,  Ypu>  the  fraction 
of  time  when  flow  is  in  the  main  stream  direction,  for  three  streamwise 
stations  in  the  wind  tunnel  (using  the  same  set  of  data  as  Figs.  5.3  and 
5.4).  Reasonably  good  agreements  are  seen  except  very  near  the  wall. 
For  x  =  317.5  cm,  intermittency  factor  is  1  throughout  the  boundary 
layer. 

If  one  divides  the  total  number  of  valid  velocity  measurements 
in  each  position  bin  to  the  total  time  that  the  measurement  volume 
spends  in  that  bin  for  the  entire  record,  T^»  one  will  have  an  idea  of 
the  data  rate  and/or  concentration  throughout  the  boundary  layer. 
Figure  5.6  shows  the  result  of  such  a  parameter  for  four  stations  in  the 
wind  tunnel.  It  is  clearly  shown  in  Fig.  5.6  that  the  data  rates  (or 
sample  size)  for  position  bins  in  the  extremes  of  the  scan  range  are 
very  low.  One  can  increase  the  data  rates  by  increasing  the  baby  powder 
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Table  5.1  Legend  for  figures  5.1  to  5.6 


x-Location 

cm. (in. ) 

Scan  frequency 

f  -  Hz. 

Counter  ampl ifier 

gain-  dB 

+Ttb 

sec. 

317.5  (125) 

55.7 

2.3  • 

2.26 

355.6  (140) 

56.6 

3.3  ♦ 

2.76 

396.2  (156) 

58.6 

3.5  A 

1.07 

436.9  (172) 

58.2 

3.9  ■ 

1 .40 

Corresponding 

sample  size  of 

pointwise  measurements 

500  each. 

using  the  same  scanning 

O 

set-up  with 

+Ttb  is  the  total  time  spent  by  the  measurement  volume  in  a  bin. 
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Y,cm. 


Fig.  5.5  Intermi ttency  Factor  vs.  Position  at  Four 
Different  Streamwise  Location. 

For  legend  see  Table  5.1  . 


Fig.  5.6  Total  Number  of  Valid  Measurements  /  Total  Time  Spent 
by  Measurement  Volume  in  a  Bin  vs.  Position. 

For  legend  see  Table  5.1  . 


94 


concentration  and/or  increasing  the  gain  of  the  counter  amplifier. 
However,  there  is  an  optimum  state  that  if  one  increases  the  gain  it 
will  amplify  the  noise.  This  can  practically  be  detected  by  the  depar¬ 
ture  of  the  mean  velocity  profile  found  by  scanning  from  the  'true' 
profile.  On  the  other  hand,  if  one  increases  the  baby  powder  concen¬ 
tration  above  certain  level,  it  either  saturates  the  PM-tube  and/or 
increases  the  noise.  In  all  of  the  experiments  reported  here  in  the 
scanning  mode,  the  following  fixed  settings  are  used: 

Laser  Power  :  600  mw 

PM-tube  Voltage  :  -1690  volts 

Average  baby 
powder  concen- 

3 

tration  :  33  mg/m 

The  scanning  system  is  very  sensitive  to  the  gain  of  the  counter 
processor  amplifier,  and  for  each  experiment  its  value  in  dB  is  written 
on  the  Table  5.1. 

As  shown  in  Fig.  3.12  section  3.2.3,  the  scanner  position 
waveform  remains  linear  up  to  156  Hz  (there  was  no  attempt  to  examine 
its  linearity  for  higher  than  156  Hz.),  and  there  is  no  mechanical 
problem.  However,  as  one  increases  the  scan  frequency,  one  should 
ensure  that  there  are  enough  number  of  fringes  crossed  by  the  particles. 
If  using  the  processing  technique  explained  here,  one  needs  faster  A/D 
and  computer,  see  section  4.7. 
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5.2  Some  Characteristics  of  the  Flow  Field  in  this  Work 

Table  5.2  shows  the  boundary  layer  integral  parameter  at  the 
same  four  streamwise  locations  using  pointwise  measurements  shown  in 
Figs.  5.3.  Flow  is  completely  attached  at  x  =  317.5  cm,  and  intermit- 
tency  factor,  Ypu,  is  1  throughout  the  boundary  layer. 

To  show  that  the  flow  is  well  behaved  in  this  work,  results 
taken  at  x  =  317.5  cm  are  compared  with  the  mean  velocity  correlation 
suggested  by  Perry  and  Schofield  (1973)  for  strong  adverse  pressure 
gradient  turbulent  boundary  layers.  Figure  5.7  shows  the  results  of  the 
non-dimensional ized  velocity  versus  normalized  position  in  which  solid 
line  is  the  correlation  suggested  by  Schofield  (1980) 

U  -  0  , 

-- -  =  1-0.4  (y/A)4  -  0.6  sin  [(*/2)  (y/A)]  (5.1) 

A 

where  Us  is  a  scaling  velocity  and  A  is  an  integral  layer  thickness 

[  =  2.86  6.  (U  /U  )].  At  this  station  U  /U  was  calculated  to  be  .99. 

L  1  '  oo'  s  s  00 

The  spectrum  function  F(n)  of  u7  is  defined  as 


(5.2) 
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Table  5.2  Boundary  layer  integral  parameters 


x-cm  (in) 

^  99~cm* 

6^-cm. 

62-cm 

H 

317.5  (125) 

10 

3.99 

1.67 

2.39 

355.6  (140) 

15 

8.34 

2.30 

3.63 

396.2  (156) 

24.5 

14.96 

3.11 

4.81 

436.9  (172) 

32.5 

19.23 

2.99 

6.43 

99 


spectra  shown  in  Fig.  5.8  are  smooth  lines  drawn  through  the  experi¬ 
mental  curves.  The  frequency  corresponding  to  the  maximum  of  the  curves 
in  Fig.  5.8  (nmax)  shifts  towards  the  lower  frequencies  as  flow  enters 
the  separation  zone.  This  means  that  the  large-scale  energy-containing 
eddies  grow  in  size  in  separated  region.  In  addition,  the  average 
frequency  characteri sties  of  the  large  energy-containing  eddies  in  the 
separated  zone  is  seen  to  be  in  the  order  of  15  Hz.  Hence,  on  the 
average,  with  60  Hz  scan  frequency,  one  has  four  scans  through  these 
eddies  which  may  reveal  more  information  about  their  structure. 

5.3 _ Quadrant  Analysis  and  Space-Time  Correlations 

The  qualitative  investigations  by  Simpson  et  al .  (1981)  showed 
that  for  the  case  where  the  thickness  of  the  backflow  region  is  small 
compared  with  the  turbulent  shear-layer  thickness;  the  backflow  region 
is  supplied  locally  by  outer  region  large-scale  structures.  The  scan¬ 
ning  LDA  system  used  here  can  be  useful  for  this  kind  of  investigation 
provided  "adequate"  data  rates  are  achievable.  The  quadrant  analysis 
and  space-time  correlation  measurements  are  attempts  toward  this  goal 
using  scanning  LDA. 

Since  the  maximum  number  of  validated  measurements,  for  a  given 
scan,  in  this  work  was  not  more  than  20,  one  needs  to  break  the  scan 
range  into  few  intervals,  Fig.  5.9,  to  be  able  to  perform  any  quadrant 
and  space-time  correlation  analysis. 

A  point  in  a  quadrant  plane  shows  the  departure  of  the  instan¬ 
taneous  streamwise  velocity  from  the  corresponding  mean  value  in  a  given 
interval,  Fig.  5.9,  versus  the  same  departure  for  another  interval  in  a 
single  scan  of  the  measurement  volume  (D-U).  Figures  5.10  to  5.16  show 


the  results  of  quadrant  analysis  using  all  scans  in  a  data  record  for 
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the  streamwise  location  of  396.2  cm.  The  number  of  points  in  each 
quadrant  is  written  in  the  corresponding  quadrant,  and  the  normalized 
space-time  correlation  is  written  in  the  parenthesis  in  the  first 
quadrant.  The  logic  in  the  quadrant  analysis  program  examines,  first, 
for  any  valid  velocity  measurements  in  the  interval  on  the  abscissa  of 
Figs.  5.10  to  5.15.  If  there  is  any,  then  the  program  searches  for 
valid  measurements  on  the  same  scan  in  the  interval  on  the  ordinate,  and 
only  if  there  is  at  least  one,  then  points  are  marked  on  the  quadrant 
plane  and  are  also  used  for  correlation  calculations  having  zero  cycle 
delay.  The  same  logic  is  used  for  measured  velocities  in  two  intervals 
having  a  pre-chosen  number  of  cycle  forward  delay,  see  Fig.  5.17,  in 
space-time  correlation  calculations. 

As  it  is  shown  in  Fig.  5.17,  one  can  find  space-time  correlation 
coefficient  between  measured  velocities  in  intervals  1  and  i  for  the 
same  cycle  throughout  the  entire  record  (zero  cycle  delay),  or  corre¬ 
lation  coefficient  between  velocities  in  the  aforementioned  intervals 
having  any  desired  number  of  forward-cycle  delay.  Theoretically, 
space-time  correlation  for  two  positions  having  t  =  time  delay  is 
defined  as: 


to41 


*WyryrT)  = 


Ju2(y1 , t)  '  u2(yi ,t) 


-  li 
T-m 


1"  /“ 


(yr  t * )  .  u(y.j ,  t '  +t  ) dt 1 


(5.4) 
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Fig.  5.10  Quadrant  Analysis  (Section  5.3). 

For  intervals  information  see  Fig.. 5. 9  . 


INTERVAL  2 


Fig.  5.11  Quadrant  Analysis(Section  5.3). 

For  intervals  information  see  Fig.  5.9 


INTERVAL  4 


Fig.  5.13  Quadrant  Ana1y$is(Section  5.3). 

For  intervals  information  see  Fig.  5.9 


Fig.  5.16  Quadrant  Analysis(Section  5.3)  . 

For  intervals  information  see  Fig.  5.9  . 
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TWO  CYCLES 


Fig.  5.17  Terminology  Used  For  Space-Time  Correlation 
Coefficients  Determined  by  the  Scanning  LDA. 


no 


In  this  work,  R^u  (y.-yj,x)  is  estimated  (Fig.  5.17)  according  to  the 
following  equation: 


Ruu<,ryi,,)  * 


J^utyj+AYj^.t)  •  u(y.±AY./2,t+T) 
^]u2(y1±AY1/2)t)  •2u2(yi±AYi/2,t+x)  * 


(5.5) 


In  which  x 


for  m  cycle  of  forward  delay,  see  Fig.  5.17, 


yryi 


summation  is  taken  over  the  entire  record.  In  Eq.  5.5,  Y.j+AY.j  in  the 
argument  of  (u)  means  that  the  value  of  this  argument  is  between 
Y.  -  AYi/2  and  y.  +  aY.j/2. 

In  all  of  the  space-time  correlation  reported  here,  a  fixed 
interval  near  the  wall  was  chosen  against  various  intervals  in  outer 
region  of  the  boundary  layer  for  a  pre-chosen  number  of  forward-cycle 
delays  to  investigate  the  influence  of  these  two  regions  on  each  other 
(quadrant  analyses  are  for  zero  cycle  delay  only,  Fig.  5.17). 

The  uncertainty  of  the  space-time  correlations  was  estimated  by 
a  comparison  between  the  auto-correlations  found  using  the  scanning  LDA 
velocity  data  and  a  fixed  measurement  volume  experiment  at  y/6  s  .5 
where  maximum  data  rate  was  available,  see  Figs.  5.9  and  5.6.  As  it  is 
shown  in  Fig.  5.18,  the  departure  of  the  scanning  results  from  that  of 
the  fixed  measurement  volume  data  increases  as  the  number  of  forward- 
cycle  delay  and/or  interval  thickness  is  increased.  For  this  reason, 
only  results  up  to  four  cycle  forward  delayed  is  reported  here.  The 
uncertainty  on  space-time  correlations  for  scanning  results  was  esti¬ 
mated  to  be  +0.1  from  the  values  calculated  using  Eq.  (5.5). 


AUTO-CORRELATION  COEFFICIENT 
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Fig.  5.18  Auto-Correlation  Found  Using  Scanning  Setup  but 
for  Fixed  Measurement  Volume  at  Y  =  1 0 . 2  cm  as 
Compared  with  the  Scanning  Results  for  Various 
Number  of  Forward  Cycle  Delay  for  an  Interval 
Centered  around  the  same  Y-Position. 

Symbol e  Interval  Jhickness;  cm._ 

*  3.4 

■  1  .7 
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The  interpretation  of  the  quadrant  analysis  is  that  when  flow  is 
greater  than  the  mean  near  the  wall  (interval  1),  then  approximately 
75±5%  of  the  time  flow  is  greater  than  the  corresponding  mean  values  in 
the  central  region  of  the  boundary  layer  as  well  (Fig.  5.12,  5.13, 
5.14);  and  when  flow  near  the  wall  is  less  than  the  mean  (larger  back- 
flow),  then  about  60  +5i  of  the  time  flow  is  greater  than  the  mean  in 
the  outer  region.  One  can  reverse  the  argument  by  saying:  given  that 
the  flow  is  greater  than  the  mean  in  the  outer  region  (conditional 
probability),  then  64  ±5%  of  the  time  flow  is  larger  than  the  mean  near 
the  wall  using  Figs.  5.12  to  5.14. 

Space-time  correlation  coefficients  for  up  to  four  forward-cycle 
delay  are  reported  here.  Figure  5.19  shows  the  correlation  coefficient 
for  zero  cycle  delay  at  three  streamwise  locations  in  the  SMU  wind 
tunnel.  As  the  flow  goes  towards  the  separation  zone  the  length  scale 
(area  under  the  curves  in  Fig.  5.19)  attains  values  of  3.55  +.88  cm,  7.9 
+  1.30  cm,  8.18  t  2.18  cm.  This  means  that  the  mean  size  of  the  eddy 
grows  downstream  which  is  in  agreement  with  the  earlier  results  from 
spectra  in  Fig.  5.8  considering  that  the  nondimensional  parameter 
U<c/6(nmax)  is  approximately  constant  (s4.71)  for  these  streamwise 
locations.  Note  that  the  ratio  of  the  length  scale  to  corresponding 
boundary  layer  thickness  is  about  .35  +  .09. 

Since  there  are  no  correlation  coefficient  measurements  avail¬ 
able  in  the  separated  region,  it  is  not  possible  to  critically  evaluate 
the  results  obtained  in  this  work.  Correlation  coefficient  measurements 
with  zero  time  delay  are 'available,  as  reported  in  a  work  by  Schubauer 
and  Klebanoff  (1951),  at  9.7  cm  upstream  of  the  separation  point  using 
hot  wire  anemometry  technique.  Figure  5.20  is  the  results  of  space-time 
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Fig.  5.19  Space-Time  Correlation  Coefficients  for  Zero  Cycle 
Delay  vs.  Position  at  Three  Different  Locations. 


10. 0 


Fig.  5.20  Space-Time  Correlation  Coefficients  between  a 
Given  Fixed  Interval  Near  the  Wall  and  Various 
Intervals  vs.  Y/5  for  Three  Streamwise  locations 

Arrows  show  the  first  interval  near  the  wall. 
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correlation  coefficients  for  three  different  streamwise  locations  and  up 
to  four  forward-cycle  delay.  The  time  delay  between  the  two  adjacent 
points  for  the  zero-cycle-delay  case  is  about  3  msec,  (use  t  given  in 
Ea.  5.5). 

The  space-time  correlations  for  y/6>.6  are  more  uncertain 
because  of  the  following  reasons.  Firstly,  the  data  rate  is  low  as  it 
is  shown  in  Fig.  5.6,  and  consequently  there  are  not  sufficient  points 
(see  Fig.  5.16)  for  proper  convergence,  and  the  total  record-time  is 
short,  also.  Secondly,  due  to  the  variation  of  the  turbulent-nontur- 
bulent  interface  and  the  fact  that  the  baby  powder  particles  tend  to 
stay  in  the  boundary  layer,  the  values  of  correlation  coefficients  are 
high  there.  Finally,  one  should  note  that  the  correlation  coefficients 
are  between  two  intervals  of  finite  length.  Fig.  5.18,  rather  than  two 
points  in  the  boundary  layer.  Note  that  this  last  problem  exists  for 
the  entire  space-time  correlation  calculations  in  this  work,  (see 
section  5.4). 

Despite  that,  one  can  see  that  the  correlation  between  the  near 
wall  region  and  the  outer  region  decreases  as  the  number  of  forward 
cycle  delay  is  increased.  It  seems  that  this  correlation  remains  longer 
in  time  for  the  streamwise  location  of  396.2  cm,  which  is  completely 
buried  in  the  separation  zone,  which  again  is  consistent  with  the  fact 
that  the  energy  containing  eddies  grow  in  size  downstream. 

Figures  5.21  to  5.23  show  the  same  information  as  given  in  Fig. 
5.20  but  in  different  format.  The  points  discussed  above  is  as  much 
physical  insight  as  possible  at  present  with  the  existing  facilities. 

Figures  5.24  to  5.28  show  the  instantaneous  velocity  profiles  as 
obtained  using  scanning  LDA  system.  It  can  be  seen  that  most  of  the 
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Fig.  5.21  Space-  Time  Correlation  Coefficients  vs.  Y/5 
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valid  measurements  are  in  the  central  third  of  the  boundary  layer. 
There  exist  scan  cycles  in  which  either  there  are  none  or  few  numbers  of 
valid  data.  Therefore,  one  cannot  really  deduce  any  pattern  using  these 
figures.  It  does,  however,  demonstrate  the  ability  of  having  almost 
instantaneous  velocity  profiles  if  sufficiently  high  data  rates  are 
achievable. 

5.4  Suggestions  for  Future  Improvements  and  Extensions 

Improvements  on  the  system  are  divided  into  three  parts: 
transmitting  optics ,  seeding  and  signal  processing.  As  for  the  trans¬ 
mitting  optics,  one  can  obtain  better  spatial  resolution  if  the  vertical 
dimension  of  the  beams  cross-over  volume  is  decreased.  One  alternative 
is  to  decrease  the  beam  diameter,  but  as  was  discussed  in  section  3.2.1 
this  means  that  better  quality  mirrors  are  needed  (flatter  surface 
curvature)  so  that  the  two  laser  beams  align  perfectly  for  all  the 
positions  of  scanning  (Fig.  3.5  c). 

Another  improvement  is  to  drive  the  scanner  mirror  with  a 
symmetric  waveform,  preferably  sinusoidal,  to  be  able  to  acquire  val¬ 
idated  measurements  for  both  upsweep  and  downsweep  scan  as  well  as 
having  better  stability  in  mirror  oscillations.  As  it  was  mentioned 
before  and  is  obvious  from  Figs.  5.6  and  5.23  to  5.27,  one  cannot  have 
more  than  a  few  consecutive  instantaneous  velocity  profiles;  and  even 
then,  there  are  not  sufficient  measurements  in  the  inner  and  outer  third 
of  the  boundary  layer.  This  problem  may  be  overcome  by  using  stair-step 
type  waveform  which  causes  the  probe  volume  to  reside  at  several  dis¬ 
crete  y-locations  and  therefore  increases  the  chance  of  having  measure¬ 
ments  for  all  the  resident  positions  in  a  given  scan.  Note  that  one 
needs  to  adjust  the  damping  ratio  of  the  scanner  control  unit  to 
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decrease  the  overshoot  of  the  step  response.  A  sinusoidal  waveform  has 
one  more  advantage  that  it  spends  more  time  near  the  wall  and  the  outer 
region  than  the  central  zone  of  the  scan  range  (from  probability  ampli¬ 
tude  distribution).  Therefore,  the  chance  of  having  more  consecutive 
"instantaneous"  profiles  would  be  increased. 

The  seeding  technique  used  here  is  quite  adequate  for  mean  and 
rms  value  measurements  and  the  only  improvement  is  to  increase  the 
record  time  which  in  turn  means  larger  memory  area  (section  5.1).  The 
seeder  itself,  though,  needs  some  improvements  to  increase  its  stability 
and/or  controllability.  There  are  two  major  improvements  on  the  seeder: 
firstly,  feed-back  control  to  maintain  a  constant  piston  velocity  for 
stability  of  the  baby  powder  concentration;  and,  secondly,  variable 
speed  wirebrush  motor  for  controlabil  ity.  In  the  absence  of  any  other 
seeding  techniques,  the  aforementioned  improvements  on  the  seeding  are 
necessary.  More  investigations  on  the  seeding  techniques  and  seeding 
materials  are  encouraged  for  any  sophisticated  flow  analysis. 

As  far  as  data  processing  is  concerned  the  counter  type  proces¬ 
sor  was  found  to  be  quite  adequate  for  this  work  despite  being  very 
sensitive  to  the  amplifier  gain  settings  and  not  having  better  scale 
sensitivity.  An  automatic  gain  control  would  be  of  some  value  in  this 
work  as  was  also  mentioned  by  Orloff  et  al .  (1975).  As  for  the  acquisi¬ 
tion  technique,  it  is  essential  to  be  able  to  sample  both  position  and 
velocity  information  simultaneously  when  the  latter  is  available  (can  be 
found  by  data  ready  pulse  from  the  counter).  The  data  acquisition  used 
in  this  work  was  the  only  choice  based  on  the  facil  ities  available  at 


Fig.  5.29  A  Suggested  Extention  of  the  Existing 

Scanning  LDA  to  Two  Dimensional  scanning  system. 
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A  transient  recorder  may  be  the  "ideal"  data  processing  tech¬ 
nique,  Durst  et  al.  (1981),  which  digitizes  Doppler  signal  from  the 
PM-tube  and  the  position  transducer  voltage  simultaneously,  and  is 
triggered  only  when  a  Doppler  burst  is  available  (Bates,  1979).  Trans¬ 
ient  recorders  with  sampling  frequency  up  to  100  MHz,  suitable  for  air 
studies,  are  available  (see  Petersen  and  Maurer,  1975). 

One  important  advantage  of  the  existing  design  is  the  possi¬ 
bility  of  extension  so  that  it  could  measure  streamwise  and  cross-flow 
velocity  profiles  simultaneously.  Figure  5.29  shows  the  schematic  of 
the  design  proposal.  The  principle  of  operation  for  this  three  beam 
system  for  point-by-point  measurement  is  discussed  in  a  report  by  TSI, 
Inc.  (1983).  The  two  side  beams  need  to  be  down  and  up-shifted  with  two 
different  frequencies,  and  path  equalization  for  three  beams  should  also 
be  considered. 

5.5  Conclusion^ 

A  rapidly  scanning  di rectional ly-sensitive  fringe-type  laser 
Doppler  Anemometer  (SLDA)  which  scans  the  measurement  volume  perpendicu¬ 
lar  to  the  optical  axis  of  the  transmitting  optics  is  described.  The 
fringe  spacing  can  be  easily  adjusted  from  1  urn  to  4^m.  Scan  fre¬ 
quencies  up  to  60  Hz  over  scan  distan-es  of  40  cm  have  been  used, 
although  scan  frequencies  up  to  150  Hz  are  possible.  The  maximum 
scanning  velocity  of  the  measurement  volume  that  can  be  used  is  directly 
proportional  to  the  shift  frequency  of  the  Bragg  cell,  since  each  signal 
producing  particle  must  cross  a  minimum  number  of  fringes  to  produce  a 
valid  signal.  Signal  averaging  oias  is  lower  with  a  scanning  IDA  than 
with  pointwise  measurements.  Results  obtained  for  a  separating  tur¬ 
bulent  boundary  layer  show  that  good  mean  and  r ms  velocity  profiles  and 
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profiles  of  the  fraction  of  time  that  the  flow  moves  downstream  can  be 
obtained  in  less  than  one  minute  of  data  acquisition.  Improvements  in 
seeding  and  signal  processing  are  needed  to  further  improve  the  data 
rate  and  permit  almost  instantaneous  velocity  profiles  to  be  obtained 
for  each  scan. 

Space-time  correlations ,  although  with  uncertainty  band,  were 
taken  for  the  first  time  in  separated  turbulent  boundary  layer  and 
showed  that  the  length  scale  for  the  large  eddies  grows  in  size  towards 
the  separation  although  the  ratio  of  the  length  scale  to  the  boundary 
layer  thickness  remains  constant.  Quadrant  analysis  not  too  far  from 
the  separation  point  says  that  when  flow  is  greater  than  the  mean  in  the 
outer  region,  then  it  is  most  probably  greater  than  the  mean  near  the 
wall,  too,  and  vice  versa,  which  represents  some  degree  of  coherence 
between  these  two  regions. 

The  advantages  of  the  SLDA  reported  here  for  large  wind-tunnel  - 
flow  surveying  is  obvious  if  one  considers  that  at  each  location  two- 
dimensional  velocity  profiles  are  possible  in  an  order  of  a  few  minutes 
(for  better  convergence).  With  improvements  mentioned  above,  one  can 
obtain  almost  instantaneous  velocity  profiles  which  could  reveal  more 
details  of  the  flow  structure. 
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